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INTRODUCTION 


The  primary  scientific  goal  for  the  Thule  optical  facility  was  to  acquire  measure¬ 
ments  of  upper  thermospheric  winds  and  temperatures  during  the  dark  observing  seasons 
covered  by  this  report:  Um  months  September  to  April  and  the  seasons  1989/90  and  1990/1. 
Unfortunately,  funding  for  the  original  three  year  effort  was  terminated  at  the  end  of  the 
second  year:  October,  1991.  Secondary  goals  were  1)  to  continue  a  lower  thermospheric 
wind  and  temperature  measurement  program  with  the  Fabry  Perot  interferometer  estab¬ 
lished  in  the  previous  report  period:  2)  to  continue  the  measurement  program  of  the  AHjL 
spectrophotometer;  3)  to  field  an  all  sky  camera  system;  4)  to  test  a  near  infrared  Michelson 
interferometer;  and  5)  to  support  various  NSF,  AFGL,  and  CEDAR  sponsored  measure¬ 
ment  campaigns. 

The  data  acquired  by  the  Fabry  Perot  interferometer  has  been  used  to  extend  the 
thermospheric  database  coverage  for  Thule  from  the  initial  December,  1985,  period  to  the 
conclusion  of  the  1990/1  observing  season.  The  initial  observations  have  been  used  in  con¬ 
cert  with  the  NCAR  Thermospheric  General  Circulation  Model  to  improve  the  understand¬ 
ing  of  the  thermodynamics  of  the  polar  cap  ionosphere  [Killeen  er  a/.,  1991;  Niciejewski  et 
al,  1992;  Niciejewski  and  Killeen,  1992;  Niciejewski  et  al,  1993;  Thayer  et  al,  1994; 
Johnson  et  al.,  1994;  Niciejewski  et  al.,  1994].  Tte  airglow  facility  at  Thule  Air  Base  is  the 
high  latitude  anchor  of  a  latitudinal  chain  of  Fabry  Perot  interferometer  stations  devoted  to 
observing  the  01  (6300  A)  emission  feature.  Thus,  the  F  region  measurements  obtained  at 
Thule  Air  Base  provide  v^uable  data  in  a  global  effort  at  improving  our  understanding  of 
the  dynamic  state  of  the  upper  atmosphere. 

The  Ebert  Fastie  spectrophotometer  has  been  used  to  observe  emissions  from  vari¬ 
ous  niehtglow  features  including  NI  (5200  A),  OI  (5577  A),  OI  (63(X)  A),  Balmer  H„ 
(6563  A),  on  (7320/30  A),  and  the  OH  (8-3)  vibrational  band.  In  addition,  several  of  these 
features  are  enhanced  in  aurora,  and  thus  the  spectrophotometer  serves  as  an  auroral  moni¬ 
tor  aiding  the  interpretation  of  lower  thermospheric  temperature  measurements.  Interpreta¬ 
tion  of  the  rotational  development  of  the  hydroxyl  band  permits  the  investigation  of  tem¬ 
perature,  gravity  wave,  and  tidal  motions  at  the  mesopause.  In  addition,  lower  atmosphere 
scattering  from  air  base  Hg  vapor  lamps  was  monitored  and  served  as  an  indicator  of  sky 
transparency  conditions. 

The  second  trailer  was  modified  by  adding  an  internal  wall  to  make  two  separate 
rooms,  and  by  adding  a  second  dome.  An  all  sky  camera  was  housed  under  a  dome  in  the 
new  trailer  during  the  1988/9  observing  season.  The  sky  was  observed  through  several  fil¬ 
ters  including  OI  (5577  A),  OI  (6300  A),  and  011  (7320/30  A)  for  the  purpose  of  coordi¬ 
nating  measurements  of  sun  aligned  ar.s  with  both  the  Fabry  Perot  interferometer  and  the 
Ebert  Fastie  spectrophotometer.  A  Michelson  interferometer  was  placed  under  the  second 
dome  and  operated  for  one  winter  season  observing  mesospheric  hydroxyl  airglow  emis¬ 
sions.  The  current  layout  of  the  Thule  facility  is  shown  in  Fi^re  1 .  Tlie  original  trailer  with 
its  two  doiiMS  is  in  tte  foreground  and  the  new  trailer  is  on  the  right. 

The  optical  facility  has  been  in  operation  during  most  of  the  winter  CEDAR  cam¬ 
paigns  during  the  report  period.  Data  has  been  acquired  for  several  LTCS  and  HLPS  cam¬ 
paigns  as  well  as  for  GITCAD/GISMOS  and  other  global  efforts. 

This  final  report  will  summarize  optical  facility  activity  in  Thule,  present  results 
from  the  varitnis  exixriments,  and  describe  woric  that  is  currently  in  progress.  Included  are 
reprints  of  papers  published  in  geophysical  journals,  and  several  abstracts  of  papers  pre- 
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sented  or  to  be  presented  at  AGU  meetings.  Two  graduate  student  theses  are  based  on  data 
acquired  at  Thule  Air  Base  and  a  report  on  the  conclusion  of  their  research  is  included. 

EXPERIMENTAL  ASPECTS  OF  THE  THULE  AIR  BASE  FACILITY 

Up  to  four  instruments  have  operated  under  separate  domes  at  the  Thule  Air  Base 
airglow  facility  in  Greenland.  The  primary  instrument  is  the  Fabry  Perot  interferometer, 
and  this  is  supplemented  by  an  Ebert  Fastie  spectrophotometer,  an  all  sky  camera  system, 
and  a  Michelson  interferometer.  The  Fabry  Perot  interferometer  has  been  observing  the  OI 
(6300  k)  emission  line  profile  since  the  start  of  operation  of  the  Thule  ground  site.  The 
scan  sequence  includes  a  measurement  in  the  zenith  followed  by  observations  in  the  four 
cardinal  directions  at  elevation  angles  of  45  degrees.  The  typical  integ^on  period  is  a  few 
minutes  for  a  1(X)  Rayleigh  signal  and  is  determined  primly  by  waiting  for  the  signal  mi¬ 
nus  background  count  to  exceed  a  predetermined  minimum.  The  statistical  uncertainties  for 
neutral  winds  and  temperatures  are  on  the  order  of  IS  m/sec  and  SO  K  respectively.  In  ad¬ 
dition  to  monitoring  the  sky  signal,  a  Ne  lamp  emission  at  6330  A  is  also  observed.  This  is 
accomplished  by  moving  the  mirror  system  to  a  position  that  looks  into  the  box  containing 
the  calibration  source.  By  calculating  a  pseudo  wind  from  the  Ne  emission,  any  drifts  in  the 
stability  of  the  instrument  may  be  h^en  into  account.  Typically,  in  any  24  hour  period,  the 
instrument  remains  stable  to  within  2S  m/sec.  Finally,  the  dark  count  rate  is  also  monitored 
and  removed  from  the  sky  signal. 

The  most  important  new  contribution  from  the  operation  of  the  Fabry  Perot  inter¬ 
ferometer  was  the  monitoring  of  the  lower  thermosphere.  This  was  accomplished  as  a  re¬ 
sponse  to  the  desire  from  the  CEDAR  aeronomical  community  to  have  measurements  of  the 
line  shape  profile  of  the  OI  (5577  A)  emission  during  the  September,  1987,  LTCS  cam¬ 
paign.  The  filter  used  has  a  bandwidth  of  3  A  centered  slightly  above  the  green  line  peak 
wavelength.  Tilting  the  filter  fine  tunes  the  Fabry  Perot's  response  to  the  sky  emission.  The 
filter  remains  in  use  in  Thule  along  with  the  OI  (6300  A)  filter  necessitating  a  bi-filter 
scanning  sequence.  At  each  pointing  direction,  first  a  green  line  fringe  is  acquit  and  then 
a  red  line  fringe  before  moving  the  mirror  system  to  a  new  direction.  With  the  new  OI 
(6300  A)  filter  mentioned  above,  a  full  scan  including  the  lower  thermosphere  measure¬ 
ments  as  well  as  the  dark  and  Ne  experiments  requires  on  the  average  20  minutes. 

The  0.5m  Ebert  Fastie  spectrophotometer  has  been  operating  at  the  Thule  site  dur¬ 
ing  the  period  of  this  report.  Prior  to  the  summer  of  1989,  the  instrument  was  interfaced  to 
the  same  data  acquisition  computer  as  the  Fabry  Perot  interferometer  and  used  the  same 
data  control  crate.  The  instrument  is  pointed  at  the  zenith  which,  in  Thule,  is  approximately 
the  geomagnetic  zenith  direction.  The  detector  is  a  GaAs  photomultiplier  tube  which  has  a 
relatively  high  quantum  efficiency  in  the  first  order  wavelengA  coverage  of  the  spectropho¬ 
tometer  -  ~38(X)  to  ~76(X)  A.  The  instrument  routinely  monitors  several  nightglow  eims- 
sions  as  well  as  the  scattering  from  low  altitude  clouds  or  fog  from  base  lifting.  During 
attended  operation,  wavelen^  and  intensity  calibrations  are  routinely  performed.  A  set  of 
spectral  lamps  permits  wavelength  calibration  while  a  calibrated  intensity  source  allows  in¬ 
tensity  calibration.  This  source  has  been  characterized  many  times  during  its  lifetime,  most 
recently  in  September,  1987.  The  source  has  a  peak  response  near  6300  A  of  7  R/A  and 
may  be  used  to  intensity  calibrate  the  spectrophotometer  over  its  full  range  of  sensitivity.  In 
addition,  while  an  operator  is  present,  the  intensity  calibration  is  transferred  to  the  Fabry 
Perot  interferometer  by  pointing  the  latter  instrument  at  the  zenith.  Knowledge  of  the  filter 
transmission  functions  plus  the  results  of  the  cross  calibration  experiments  allows  a  reason¬ 
able  estimate  as  to  the  true  brightness  of  the  nightglow  viewed  by  the  Fabry  Perot  interfer¬ 
ometer. 

In  an  effort  to  increase  the  efficiency  of  tlK  the  spectrophotometer  experiment,  the 
instrument  was  interfaced  with  an  IBM  AT  clone  computer  during  the  1989  down  season. 
This  was  accomplished  at  SPRL  by  two  students  funded  by  the  REU  program.  The  ra- 
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tionale  for  this  major  change  was  twofold:  one.  to  make  the  experiments  in  Thule  truly  in¬ 
dependent;  and  two,  to  upgrade  the  computer  control  system  to  current  state  of  the  an 
levels.  The  upgrade  was  motivated  by  the  successes  that  were  made  during  the  startup  ef¬ 
forts  in  September,  1987,  and  November,  1988.  Both  the  Fabiy  Perot  interferometer  and 
the  Ebert  Fastie  spectrophotometer  experiments  performed  so  well  that  the  single  data  ac¬ 
quisition  computer  became  overloaded  with  data  prior  to  a  site  visit  and,  as  a  result,  shut 
down  for  lack  of  disk  space.  Rather  than  increasing  disk  space  on  the  old  PDP  1 1/23  com¬ 
puter,  weaning  the  spectrophotometer  experiment  away  from  this  computer  and  interfacing 
with  its  own  system  with  an  optical  disk  drive  data  storage  unit  was  the  preferred  solution. 
This  new  arrangement  was  tested  in  Ann  Arbor  in  September,  1989,  and  was  shipped  and 
installed  at  Thule  Air  Base  in  November,  1989. 

During  the  fall  of  1987,  an  additional  trailer  was  placed  next  to  the  original  Wells 
Fargo  trailer.  The  new  trailer  initially  had  one  hole  cut  into  the  ceiling  with  a  clear  Plexi¬ 
glass  dome  placed  over  it.  In  November,  1988,  an  all  sky  camera  system  was  placed  under 
this  dome.  The  detection  system  for  this  experiment  makes  use  of  a  bare  CCD  placed  at  the 
focal  plane  of  an  all  sky  camera  lens  system.  The  filter  wheel  with  its  complement  of  filters 
was  placed  above  the  input  optics.  Since  it  was  decided  at  the  very  start  that  there  would 
not  be  any  image  intensification,  the  CCD  must  be  kept  as  cold  as  possible  to  reduce  the 
dark  count  inherent  within  the  device.  This  is  accomplished  in  several  ways.  First  of  all, 
the  trailer  housing  this  experiment  is  kept  at  ambient  tenq)erature  which  is  near  -30**  C  dur¬ 
ing  the  winter  months.  Secondly,  the  Peltier  cooling  system  that  was  supplied  by  the  ven¬ 
dor  was  supplemented  with  a  thermoelectric  cooling  system.  Dark  counts  were  observed  to 
be  as  low  as  S  counts  per  300  seconds.  The  current  version  of  the  experiment  is  capable  of 
imaging  sub-visual  auroral  arcs  and  patches  at  5577  and  6300  A.  The  images  are  recorded 
on  video  tape  for  playback  and  analysis  in  Michigan. 

During  the  fall  of  1990,  a  Michelson  interferometer  was  placed  under  a  new  second 
dome  in  the  second  trailer.  The  instrument  is  controlled  by  a  separate  computer  and  runs 
independently  of  the  three  other  experiments,  pie  bandwidth  observe^  by  the  Michelson 
interferometer  is  the  near  infrared  from  lOOOOA  to  I6OOOA  which  includes  many  intense 
hydroxyl  airglow  emissions  as  well  as  two  bands  of  the  near  infrared  O2  IR  Atmospheric 
band.  The  instrument  was  setup  to  obtain  two  minute  integration  spectra  during  twilight 
and  ten  minute  integrations  during  dark  conditions.  Figure  2  shows  an  example  of  a  re¬ 
duced  data  set  for  one  night  at  Thule  Air  Base.  The  top  half  shows  a  time  series  of  me- 
sopause  neutral  temperatures  extracted  from  the  Meinel  OH  (3,1)  band.  The  bottom  half 
shows  the  corresponding  time  series  of  intensity  variations. 
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RESULTS 


The  Thule  Air  Base  facility  was  brought  on  line  in  November,  1989  by  Dr.  Ni- 
ciejewski  and  graduate  students  Matt  Turnbull  and  Elaine  Tnidell.  Dr.  Niciejewski  initiated 
the  FPI  and  led  the  effort,  Mr.  Turnbull  initiated  the  EPS,  while  Ms.  Trudell  initiated  the 
all-sky  camera  system  operation.  The  site  was  again  visited  by  Dr.  Niciejewski  in  Febru¬ 
ary,  19^  to  back  up  data  sets  as  well  as  to  fix  a  software  problem  with  the  EFS.  Mr. 
Turnbull  was  charged  with  shutting  down  the  station  for  the  summer  term  as  well  as  return¬ 
ing  the  ^1-sky  camera  to  Michigan  during  a  visit  in  April,  1990.  The  site  was  next  visited 
by  Andy  Betki  during  August,  1990  to  perform  modiEcations  to  the  second  trailer.  The 
1990/1  observing  season  was  initiated  by  a  visit  in  September,  1990,  by  Mr.  Turnbull  and 
Drs.  Niciejewski  and  Murty.  During  that  trip,  a  Michelson  interferometer  was  added  to  the 
instrumental  complement.  Unfortunately,  the  photomultiplier  tube  for  the  EFS  did  not 
function  so  the  EFS  was  not  brought  on  line  at  t^t  time.  The  site  was  visited  by  Dr.  Murty 
in  December,  1990  to  back  up  old  FPI  and  MI  data.  The  site  was  next  visited  by  Dr.  Ni¬ 
ciejewski,  Capt.  Stephen  Carr,  and  Ms.  Trudell  in  January,  1991.  The  EFS  was  repaired 
and  placed  in  operation  for  the  rest  of  the  season.  Ms.  Trudell  was  charged  with  installing 
the  ^1-sky  camera  experiment.  Unfortunately,  one  box  of  the  shipment  did  not  arrive  on 
time  which  meant  that  this  experiment  did  not  operate  during  that  season.  Dr.  Killeen  vis¬ 
ited  the  site  in  February,  1991,  to  back  up  old  data  sets.  Shutdown  was  performed  in  April, 
1991,  by  Capt.  Carr,  who  also  returned  the  MI  and  the  all-sky  camera  boxes  to  Michigan. 

During  the  periods  mentioned  above,  high  quality  data  were  acquired  by  the  FPI 
between  January,  1990  and  February,  1990;  January,  1991  and  February,  1991.  The  EFS 
acquired  data  between  February,  1990  and  March,  1990;  January,  1991  and  March,  1991. 
The  MI  acquired  data  between  January,  1991  and  March,  1991 .  Cloud  cover  records  have 
been  acquired  from  the  Thule  Air  Base  meteorological  station.  Much  of  the  Fabry  Perot  in¬ 
terferometer  data  and  all  the  cloud  cover  records  have  been  entered  into  an  "in  house"  data 
base  at  SPRL.  Figure  3  shows  some  results  obtained  from  the  data  base.  The  average  hori¬ 
zontal  neutral  wind  meridional  and  zonal  components  at  F  region  altitudes  are  shown  for 
conditions  where  Kp  ^  3  (top)  and  3  <  Kp  <  6  (bottom).  Figure  4  displays  averaged  com¬ 
ponents  for  E  region  altitudes. 

Several  papers  resulting  from  the  analysis  of  Thule  Fabry  Perot  interferometer  data 
have  been  published  during  this  report  interval.  The  first  paper  [Killeen  etai,  1991]  was 
the  result  of  a  contribution  to  an  lUGG  conference  and  described  Thule  Air  Base  measure¬ 
ments  in  the  context  of  other  high  latitude  measurements.  The  second  paper  [Niciejewski  et 
al,  1992]  described  the  effect  of  the  orientation  of  the  By  component  of  the  IMF  on  neutral 
winds  at  high  latitudes.  The  third  paper  [Niciejewsld  and  Killeen,  1992]  described  the  ef¬ 
fect  of  Kp  on  the  magnitude  of  the  high  latitude  meridional  neutral  wind  flow  and  described 
the  appearance  of  midnight  abatement  features  in  response  to  extremely  low  Kp  conditions. 
The  fourth  paper  [Niciejewsld  et  al,  1993]  described  observations  of  neutral  winds  in  the 
presence  of  sun-aligned  auroral  arcs  above  Thule  Air  Base.  This  paper  is  currently  under¬ 
going  revision.  The  fifth  paper  [Thayer  et  al,  1994]  described  an  experimental  ground- 
based  measurement  of  the  ion/neutrd  coupling  at  very  high  geomagnetic  latitudes.  The 
sixth  paper  [Johnson  et  al.,  1994]  described  observations  and  model  predictions  of  neutral 
winds  during  the  January,  1988  GISMOS  experimental  period.  The  seventh  paper 
[Niciejewski  et  al,  1994]  described  the  effect  of  the  orientation  of  the  B^  component  of  the 
IMF  on  neutral  winds  at  high  latitudes. 

Many  presentations  at  scientific  conferences  have  been  made  during  the  report  pe¬ 
riod  which  include  data  sets  obtained  at  the  Thule  Air  Base  optical  facility.  One  paper  was 
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presented  at  the  Fall  '89  AGU  meeting.  Three  papers  were  presented  at  the  Spring  '90 
AGU  meeting.  One  paper  was  presented  at  the  COSPAR  '90  meeting.  One  paper  was  pre¬ 
sented  at  the  Spring  '91  AGU  meeting.  One  paper  was  presented  at  Ae  URSI  '91  meeting. 
Three  papers  were  presented  at  the  lUGG  '91  meeting.  One  paper  was  presented  at  the  Fall 
'91  AGU  meeting.  Two  papers  were  presented  at  the  Spring  '92  AGU  meeting.  One  paper 
was  presented  at  the  COSPAR  '92  meeting.  Three  papers  were  presented  at  tte  1992  AGU 
Chapman  conference.  Two  papers  were  presented  at  the  Spring  '93  AGU  conference. 

The  facility  at  Thule  Air  Base  has  been  used  to  train  graduate  students  in  the  opera¬ 
tion  and  maintenance  of  ground  based  optical  instmmentation.  During  the  period  of  this  re¬ 
port,  the  station  has  been  visited  by  Mr.  M.  Turnbull,  Capt.  S.  Carr,  and  Ms.  E.  Trudell. 
In  addition,  Mr.  D.  Drob  and  Mr.  Y.  Won  have  been  involved  in  the  data  analysis  and  in¬ 
terpretation  of  Thule  data.  Mr.  J.  Thayer,  Mr.  Y.  Won,  and  Capt.  S.  Carr  have  completed 
PhD  degrees  using  data  sets  accumulated  at  Thule  Air  Base.  Mr.  Drob  and  Ms.  Tmdell  are 
completing  their  graduate  degree  requirements  using  data  obtained  from  the  Thule  site,  the 
forn^r  using  MI  and  FPI  data,  the  latter  using  FPI  and  all-sky  camera  data. 

SUMMARY 

A  computer  controlled,  automatic  airglow  facility  has  been  in  operation  at  Thule  Air 
Base  during  the  winter  observing  periods  covered  by  this  report.  The  instrumentation 
complement  includes  a  Fabry  Perot  interferometer  monitoring  thermospheric  emissions 
from  OI  (5577  A)  and  01  (6300  A),  a  spectrophotometer  monitoring  airglow  emissions,  an 
all  sky  imaging  system  monitoring  visual  to  sub-visual  auroral  arcs  and  patches,  and  a  near 
infrared  Fourier  transform  spectrometer  monitoring  hydroxyl  emissions.  Thermospheric 
neutral  winds  and  temperatures  have  been  measured  routinely  since  1985  leading  to  an  ini¬ 
tial  description  of  polar  cap  thermodynamics  during  the  solar  minimum  and  solar  maximum 
phases. 
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FIGURE  CAPTIONS 


Figure  1.  The  Thule  optical  facility  consists  of  two  instrumented  trailers  for  scientific  re¬ 
search.  The  trailer  in  the  foreground  on  the  left  houses  the  Fabry  Perot  interferometer  and 
Ebert  Fastie  spectrophotometer.  The  trailer  on  the  right  houses  the  all  sky  camera  system 
and  the  Michelson  interferometer,  the  former  under  a  small  dome  not  seen  in  the  photo¬ 
graph. 

Figure  2.  Summary  plot  of  mesopause  data  for  the  date  09/24/90  as  acquired  by  the  Thule 
Air  Base  Michelson  interferometer  experiment.  The  top  figure  displays  a  time  series  of 
neutral  temperatures  as  derived  from  a  series  of  Meinel  OH  (3,1)  band  spectra.  The  bottom 
figure  displays  the  corresponding  OH  intensity  time  series. 

Figure  3.  Component  plot  of  averaged  meridional  and  zonal  neutral  winds  in  the  upper 
thermosphere  as  observed  over  Thule  Air  Base  for  conditions  of  Kp  ^  3  (top)  and  3  <  Kp  < 
6  (bottom). 

Figure  4.  Component  plot  of  averaged  meridional  and  zonal  neutral  winds  in  the  lower 
thermosphere  as  observed  over  Thule  Air  Base  for  the  winter  of  1988/89  (top)  and  the 
winter  of  19%/91  (bottom). 
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Abstract — Recent  experimental  measurements  of  the  dynamics  of  the  neutral  upper  thermosphere  have 
demonstrated  the  important  roles  of  ion-drag  and  Joule  heating  processes  in  establishing  the  basic  neutral 
wind  morphology  and  controlling  neutral  composition,  particularly  in  the  high-latitude  region.  Instruments 
on  the  Dynamics  Exptorer-2  spacecraft  (DE  2),  for  example,  were  capable  of  measuring  the  three-dimen¬ 
sional  vector  tMutral  wind  and  ion  drift  in  the  thermosphere  along  the  orbital  track,  together  with 
constituent  densities  and  temperatures.  Ground-based  optical  and  radar  measurements  of  winds  and 
temperatures  from  observatories  in  Greenland  have  contributed  additional  measurements  of  thermospheric 
neutral  wind  velocities  and  ionospheric  parameters.  The  comprehensive  nature  of  these  various  data  sets 
has  enabled  more  stringent  experimental  constraints  to  be  placed  on  the  numerical  models  of  the  region 
(thermosphere-ionosphere  general  circulation  models.  TiGCMs).  leading  to  an  improved  theoretical 
understanding  of  the  important  physical  processes  that  control  thermospheric  circulation  and  variability. 
In  addition,  the  measurements  have  enaUed  the  development  of  semi-empirkal  models  of  thermosphere 
dynamics  which  can  be  used  in  various  theoretical  studies.  The  Vector  Spherical  Harmonic  (VSH)  model, 
for  example,  provides  a  description  of  global  thermospheric  state  variables  (wind,  temperature  and  density), 
using  a  combination  of  empirical  data  and  NCAR-TIGCM  calculations.  This  paper  presents  a  brief  review 
of  some  of  the  more  recent  progress  made  in  this  area  by  the  team  of  researchers  at  the  University  of 
Michigan,  with  emphasis  on  the  interpretation  of  experimental  measurements  made  from  DE2  and  from 
ground-based  observatories  in  Thule  and  Sondrestromfjord.  Greenland.  Comparisons  between  individual 
data  sets  from  these  sources  and  the  VSH  model  arc  also  presented. 


I.  INTRODUCTION 

The  dynamics,  thermodynamics,  and  compositional 
structure  of  the  neutral  gas  in  the  upper  thermosphere 
are  strongly  controlled  by  ion-neutral  collision  pro¬ 
cesses.  particularly  at  high  geomagnetic  latitudes. 
Magnctospheric  convection  electric  fields  map  down 
into  the  high-latitude  ionosphere  along  equipotential 
geomagnetic  field  lines  and  drive  the  charged  particles 
there  into  cellular  motion.  Rapidly  moving  ions  in 
these  ionospheric  convection  cells  can  readily  transfer 
momentum  to  the  neutral  gas  via  ion-drag  forcing  and 
can  convert  a  portion  of  their  kinetic  energy  of  motion 
into  internal  energy  of  the  neutral  gas  via  Joule  or 
collisional  heating.  In  addition,  the  enhanced  con¬ 
centrations  of  ions  associated  with  auroral  particle 
precipitation  serve  to  increase  the  rates  of  exchange 
of  energy  and  momentum  between  the  neutral  and 
ionized  species  by  increasing  the  magnitude  of  the  ion- 
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neutral  collision  frequency  amongst  other  factors.  The 
effect  of  these  high-latitude  ion-neutral  collision  pro¬ 
cesses  is  to  modify  the  neutral  wind,  temperature  and 
compositional  structure  away  from  the  structure  that 
would  be  expected  if  the  only  sources  of  energy  for 
the  thermosphere  were  solar  UV  and  EUV  insolation. 

Many  examples  of  the  effects  of  auroral/ 
magnetospheric  forcing  of  the  neutral  thermosphere 
at  high  latitudes  have  been  reported  in  the  literature 
over  the  past  10  years  (e.g.  see  the  review  articles 
of  Meriwether,  1983;  Roble,  1983;  Mayr  ei  at., 
1985;  Killeen,  1987;  Killeen  and  Roble,  1988).  Ra¬ 
pid  experimental  and  theoretical  progress  has  been 
sustained  due.  in  part,  to  the  comprehensive  nature 
of  new  data  sets  provided  by  the  NASA  Dynamics 
Expiorer-2  spacecraft  (DE2)  and  the  ground- 
based  network  of  optical  interferometer  observatories 
and  incoherent  scatter  radar  facilities,  as  well  as  to 
the  maturity  of  the  numerical,  three-dimensional, 
time-dependent,  thermosphere/ionosphere  general 
circulation  models  (TIGCMs).  The  two  best-developed 
examples  of  the  latter  are  the  National  Center  for 
Atmospheric  Research  model  (NCAR-TIGCM)  of 
Dickinson  et  at.  (1981)  and  Roble  et  at.  (1982. 
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1988)  and  Ihc  University  College  London/Shefficid 
University  model  of  Fuller-Rowkll  and  Ria-s  (1980) 
and  FuLLER-Rowt-LL't’/ «/.  (1987).  Two  semi-empiri¬ 
cal  models  of  neutral  thermosphere  dynamics  have 
also  been  developed  in  recent  years,  the  Vector  Spheri¬ 
cal  Harmonic  (VSH)  model  or  Killeen  ei  al.  (1987) 
and  the  Horizontal  Wind  Model  (HWM)  of  Heoin  et 
al.  (1988.  1991).  While  the  HWM  model  is  entirely 
based  on  observations,  the  VSH  model  is  a  hybrid, 
based  primarily  on  the  theoretical  calculations  from 
the  NCAR-TIGCM.  but  supplemented  with  DE2 
measurements  where  available. 

In  this  paper,  we  summarize  some  of  the  recent 
high-latitude  results  obtained  using  the  DE2  ther¬ 
mospheric  neutral  wind,  temperature  and  com¬ 
position  data,  together  with  VSH  model  predictions. 
Additional  experimental  results  from  ground-based 
optical  observatories  sited  at  Thule  and  Son- 
drcstramQord  (Sondrestrom).  Greenland  ’•ni  from 
the  incoherent  scatter  radar  at  Sondrestrom  are  also 
shown.  The  emphasis  of  this  discussion  is  on  per¬ 
turbations  to  the  solar-driven  upper  thermospheric 
winds,  temperatures  and  compositional  structures  due 
to  forcing  associated  with  high-latitude  ionosphere- 
thermosphere  coupling  processes.  In  the  following 
sections  we  ( I )  discuss  the  nature  of  the  experimental 
measurements  made  on  the  DE2  spacecraft  and 
describe  some  recent  results.  (2)  describe  briefly  the 
VSH  model,  some  sample  results,  and  experiment- 
model  comparisons  using  both  satellite  and  ground- 
based  measurements  of  high-latitude  thermosphere 
dynamics.  (3)  discuss  considerations  pertaining  to 
compositional  perturbations  at  high  latitudes,  and  (4) 
make  some  concluding  remarks. 

2.  DYNAMICS  F-XPIXIRER  MF.A.SURF.MF,NTS 

The  Dynamics  Explorer-2  spacecraft  payload 
included  the  Fabry-Perot  interferometer.  FPI  (Hays 
el  al.,  1981).  the  Wind  and  Temperature 
Spectrometer.  WATS  (Spencer  et  al.,  1981).  and  the 
Neutral  Atmosphere  Composition  Spectrometer. 
NACS  (Carionan  et  al.,  1981).  These  instruments 
measured  the  meridional  and  zonal  components  of 
the  neutral  wind  and  neutral  constituent  abundances, 
respectively,  along  the  track  of  the  polar-orbiting  sat¬ 
ellite.  Both  the  FPI  and  WATS  instruments  also  mea¬ 
sured  neutral  kinetic  temperatures.  In  addition  to 
these  instruments,  the  Langmuir  probe,  LANG 
(Krehbiel  et  al.,  1981),  the  Ion  Drift  Meter,  IDM 
(Heelis  el  al.,  1981).  and  the  Retarding  Potential 
Analyzer,  RPA  (Hanson  et  al.,  1981)  measured  the 
ion  density  and  the  zonal  meridional  components  of 
the  ion  drift,  respectively.  The  RPA  and  LANG 


instruments  also  enabled  measurements  of  ion  and 
electron  temperatures,  respectively.  The  com¬ 
prehensive  nature  of  the  DE  2  data  set  has  enabled 
various  studies  of  thermospheric  ion  -neutral  coupling 
at  high  latitudes  to  be  conducted  (sec  the  review  by 
Killeen  and  Roble.  1988). 

Examples  of  the  DE  2  coverage  for  two  nearly  full 
orbits  of  the  spacecraft  are  shown  in  Figs  I  and  2. 
These  figures  show,  as  a  function  of  various  geo¬ 
physical  parameters  along  the  orbital  track  of  the 
spacecraft,  the  following  measured  or  derived  observ¬ 
ables  (from  top  to  bottom) ;  (I)  the  neutral  wind  vec¬ 
tor  from  FPI  and  WATS ;  (2)  the  ion  drift  vector  from 
IDM  and  RPA;  (3)  the  atomic  oxygen,  molecular 
nitrogen  and  electron  number  densities  from  NACS 
and  LANG,  respectively ;  (4)  the  electron,  ion.  and 
neutral  kinetic  temperatures  from  LANG.  RPA  and 
FPI,  respectively :  and  (S)  vertical  winds  and  ion  drifts 
measured  from  WATS  and  IDM.  respectively.  The 
upper  inset  to  the  left  of  each  figure  is  a  dial  showing 
the  neutral  winds  (Fig.  2)  or  ion  winds  (Fig.  I)  in  a 
geographic  polar  projection  for  the  southern  hemi¬ 
sphere.  The  lower  inset  provides  information  on  the 
time  history  of  the  relevant  interplanetary  magnetic 
field  (IMF)  components  and  the  and  AE  geo¬ 
magnetic  indices  prior  to  and  at  the  time  (given  by  the 
dotted  vertical  line)  of  the  orbital  pass.  The  vectors 
arc  plotted  such  that  local  noon  is  to  the  lop  of  the 
diagram  and  local  dawn  to  the  right.  Thus,  the  down¬ 
ward-directed  ion  drift  vectors  arc  in  the  anti-sunward 
direction  and  arc  generally  associated  with  Ihc  central 
region  of  the  northern  and  southern  hemisphere  geo¬ 
magnetic  polar  caps,  while  the  regions  of  upward- 
directed  ion  drift  vectors  arc  in  the  sunward  direction 
and  arc  generally  associated  with  the  dawn  and  dusk 
sectors  of  the  auroral  zones. 

Figure  I  presents  data  from  orbit  7366  of  DE2 
which  occurred  during  Julian  day  82338.  As  can  be 
seen  from  the  lower  inset,  this  orbit  occurred  dur¬ 
ing  relatively  quiet  geomagnetic  conditions,  with 
K^=  l-t-(/4,  =  13).  By  contrast.  Fig.  2  shows  data 
from  orbit  7219  on  day  82328  which  was  a  period  of 
intense  geomagnetic  activity,  with  Kp=7  (A^  =  83). 
In  both  cases,  the  satellite  obtained  data  across  both 
northern  (winter)  and  southern  (summer)  high-lali- 
ludc  regions,  as  well  as  low-  and  mid-latitudes.  These 
data  arc  shown  to  provide  an  overview  of  the  vari¬ 
ations  in  thermospheric  state  variables  typically 
observed  from  DE  2. 

There  are  several  general  features  that  can  be  dis¬ 
cerned  from  Figs  I  and  2.  Most  obviously,  the  neutral 
and  ion  winds  (including  the  horizontal  and  vertical 
components)  have  their  greatest  structure  and  mag¬ 
nitudes  in  the  two  polar  regions,  with  clear  evidence 
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Fig.  I.  Geophysical  observables  measuml  along  the  track  of  Dynamics  Explorer  2  during  orbit  7.766.  The 
neutral  winds  and  ion  drifts  are  shown  in  the  lop  two  traces  plotted  against  time,  altitude,  and  latitude  of 
the  spauecrafl  as  it  succesitivcly  passes  (len-tevrighi)  over  the  northern  hemisphere  (winter)  polar  region, 
the  equatorial  region  and  the  southern  hemisphere  (summer)  polar  region.  The  second  panel  shows  the 
atomic  oxygen  and  molecular  nitrogen  densities,  referred  to  a  fiducial  altitude  of  .100  km,  and  the  electron 
(ion)  density.  The  third  panel  shows  the  electron  (dotted  line),  ion  (light  line)  and  neutral  (bold  line) 
temperatures  measured  along  the  orbital  track.  The  bottom  trace  shows  the  measured  vertical  neutral 
winds  (bold  line)  and  ion  drifts  (light  line).  The  lower  inset  to  the  left  shows  the  lime  history  of  the  IMF 
components  and  the  geophysical  indices  AE  and  for  the  24  h  prior  to  the  orbital  pass.  The  mid-point 
UT  for  the  satellite  pass  is  denoted  by  the  vertical  dotted  line.  The  upper  inset  to  the  Icfl  shows  a  polar 
dial  (geographic  latitude  and  solar  local  time)  for  the  pass,  with  the  ion  winds  plotted  and  the  location  of 

the  geomagnetic  pole  given  by  the  cross. 


of  a  predominantly  two-cell  ion  convection  pattern 
and  a  morphologically  similar,  but  less  strong,  neutral 
wind  system  at  high  latitudes.  The  neutral  winds  (and 
ion  drifts)  at  mid-  and  low-latitudes  (central  region 
of  figures)  have  much  smaller  magnitudes  and  are 
consistent  with  a  day-to-night  solar-driven  global 
thermospheric  wind  system.  The  neutral  composition 
measurements  obtained  during  these  passes  have  all 
been  adjusted  to  an  altitude  of  300  km  in  order  to 
highlight  latitudinal  rather  than  altitudinal  changes, 
using  a  diffusive  equilibrium  approximation  to  map 
the  measured  values  upwards  or  downwards  from  the 
corresponding  spacecraft  altitude  for  each  measure¬ 
ment.  The  [O]  and  [N:]  densities,  referred  to  300  km 
altitude,  show  relatively  smooth  and  unstructured 
behavior  at  mid-  and  low-latitudes,  and  more  struc¬ 


tured  perturbations  at  polar  latitudes,  in  response  to  the 
temporally  and  spatially  complex  high-latitude  ion 
drag  and  Joule  heating  processes.  For  both  the  quiet¬ 
time  orbit  (7366)  and  the  active-time  orbit  (7219),  the 
N,  densities  at  300  km  equal  or  exceed  the  cor¬ 
responding  O  densities  in  the  summer  hemisphere 
polar  regions  due  to  the  extra  heating  and  consequent 
upwelling  of  the  thermosphere.  The  O  densities  do 
not  show  the  strong  summer-to-wintcr  gradient 
observed  in  the  N2  densities  for  the  active  period  and, 
in  fact,  show  decreases  at  high  latitudes  during  the 
disturbance  when  N,  increases.  These  ‘anomalies'  are 
well  understood  and  are  related  to  the  different  masses 
of  the  species.  It  is  also  interesting  to  note  that  the 
vertical  velocities  tend  to  be  more  upwards  in  the 
warmer  summer  hemisphere  and  in  regions  of 
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Fig.  2.  Geophysical  observables  measured  along  the  track  of  Dynamics  Explorer  2  during  orbit  7219.  The 
formal  of  the  figure  is  as  for  Fig.  I,  with  the  exception  that  neutral  winds  rather  than  ion  drifts  are  shown 

in  the  polar  dial  to  the  upper  left. 


enhanced  auroral  heating.  The  electron  density  (n,)  hemisphere  polar  passage  for  orbit  7219,  where  the 

measurements  show  various  structures  associated  greatly  enhanced  frictional  heating  of  the  ions,  due 

with  ionospheric  morphology.  A  particular  feature  of  to  the  large  ion-neutral  difference  velocities,  caused 
interest  occurs  near  68200  sUT,  where  a  depression  localized  enhancements  in  ion  temperature  to  many 

in  n,.  at  equatorial  latitudes  is  probably  an  altitude  thousands  of  degrees.  In  general,  all  three  species  are 

effect  caused  by  the  relatively  high  F-region  peak  and  in  thermal  equilibrium  at  low  latitudes  {T,.  %  T, «  T„), 

the  relatively  low  satellite  altitude;  in  this  region  with  major  departures  from  that  equilibrium  occur- 

the  satellite  dips  below  the  F;  peak.  The  ver-  ring  at  high  latitudes. 

tical  wind  measurements  (bottom  panel)  show  evi-  The  effects  of  enhanced  levels  of  geomagnetic 

dence  for  the  presence  of  gravity  wave  activity  at  activity  can  be  seen  by  comparing  the  data  of  Fig.  I 

polar  latitudes  with  various  horizontal  wavelengths,  and  Fig.  2.  The  neutral  and  ion  winds,  neutral  densi- 

ranging  from  a  few  degrees  to  tens  of  degrees  in  ties,  and  neutral  and  ion  temperatures  are  typically  all 

latitude.  greater  in  magnitude  during  the  more  geomagnetically 

The  measured  ion,  electron  and  neutral  tern-  active  orbit  7219  than  during  orbit  7366.  Neutral 

peraturesITr,  r,.  and  r„,  respectively)  along  the  orbital  winds  in  excess  of  1200  m/s  were  observed  near 

tracks  provide  insight  into  the  thermal  structure  of  67400  s  LIT  on  orbit  7219,  which  occurred  during  an 

the  upper  atmosphere.  As  expected,  the  magnitude  of  extended  period  of  greatly  enhanced  geomagnetic 

T,  is  the  greatest  of  the  three,  except  for  localized  activity  (see,  for  example,  the  AE  trace  in  the  lower 

regions  of  very  strong  Joule  heating  where  can  inset  of  Fig.  2).  These  are  the  largest  magnitude  winds 

exceed  T,.  Such  a  period  occurred  during  the  southern  observed  by  DE  and  it  is  our  belief  that  these  are  the 
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Rg.  3.  Smultaneously  measured  neutral  wind  vectors  from  DE2  (orbit  1813)  and  DE 1  auroral  image 
(courtesy  of  L.  A.  Frank  and  J.  D.  Craven,  University  of  Iowa).  The  images  were  obtained  using  the  Spin- 
Scan  Auroral  Imager,  SAI,  on  Dynamics  Explorer  1 .  viewing  the  northern  hemisphere  aurorae  at  ultraviolM 
waveiengths.  The  ""ap*  false-color  coded  and  ate  oriented  such  that  the  direction  towa^  the  Sun  is 
to  the  top  of  the  figure,  dusk  to  the  left.  The  solar  terminator  is  evident,  running  roughly  horizontal  across 
oar-ii  image,  as  is  the  entire  auroral  oval  located  just  to  the  ni^tside  of  the  terminator.  The  neutral  wind 
vectors  ate  denoted  by  the  yellow  arrows  whose  origins  are  positioned  along  the  DE  2  orbital  track.  The 
wind  scale  is  given  at  lower  right.  3b-d  are  for  DE  2  orlnts  1815,  1819  and  1823,  respectively.  The  figure 

is  taken  from  Killees  el  aJ.  (1988). 
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Fig.  7.  Neutral  winds  and  ion  drift  measurements  (blue  arrows  and  red  bars,  respectively)  are  plotted  along 
the  orbital  track  of  Dynamics  Explorer  2  during  orbits  7210  and  7211.  The  data  are  plotted  in  geomagnetic 
polar  coordinates  (magnetic  local  time  and  magnetic  latitude),  with  the  outer  circle  being  at  40°  geomag¬ 
netic.  These  data  were  obtained  during  a  prolonged  period  of  strongly  northward  IMF.  The  figure  is 

adapted  from  Killeen  et  a/.  (198S). 
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Fig.  10.  Thermospheric  wind  and  temperature  plots  for  December  solstice  solar  maximum  conditions,  (a) 
VSH  model  predictions:  (b)  averaged  DE3  wmda  and  lempemtures  coUected  between  November  1^1 
and  January  1982  and  between  November  1982  and  Jantiaty  1983.  The  plots  are  in  geographic  latitude 

and  Universal  Time. 
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largest  magnitude  winds  yet  reported  for  any  part  of 
the  Earth’s  atmosphere. 

While  the  study  of  thermospheric  and  ionospheric 
parameters  from  individual  high-latitude  orbital 
passes  provides  information  on  the  response  of  the 
upper  atmosphere  to  energy  and  momentum  inputs 
from  the  magnetosphere,  the  study  by  Killeen  et 
al.  (1988)  illustrated  unambiguously  the  relationship 
between  thermospheric  winds  and  one  of  the  most 
obvious  manifestations  of  magnetospheric  pro¬ 
cesses — the  visible  aurora.  Figure  3  shows  simul¬ 
taneous  measurements  of  global-scale  auroral 
luminosity  distributions  and  vector  neutral  winds  over 
the  northern  (winter)  polar  cap,  using  data  from  the 
spin-scan  auroral  imager,  SAI  (Frank  et  al.,  1981), 
on  DE 1  and  the  FPl  and  WATS  instruments  on  DE  2, 
respectively.  This  work  illustrates  the  spatial  relation¬ 
ship  between  large-scale  morphological  features  of  the 
F-region  neutral  wind  field  in  the  winter  polar  region 
and  the  location  and  spatial  extent  of  the  aurora. 
A  definite  correlation  can  be  seen  to  exist  between 
reversals  and  boundaries  in  the  neutral  wind  field  and 
the  location  of  the  visible  auroral  oval.  Examples  of 
such  simultaneous  data  sets  indicate  that  the  neutral 
wind  and  boundaries  closely  follow  the  substorm- 
dependent  expansion  and  contraction  of  the  auroral 
oval. 

2.1.  Dependence  of  Mgh-latitude  thermospheric  circula¬ 
tion  on  magnetospheric/ionospheric  forcing  processes 

To  study  the  dynamic  response  of  the  high-latitude 
thermosphere  to  various  forcings  in  greater  detail, 
data  from  many  hundreds  of  DE  2  orbits,  similar  to 
those  depicted  in  Figs  1  and  2,  have  been  analysed. 
The  FPI  and  WATS  measurements  of  the  two  com¬ 
ponents  of  the  upper  thermospheric  horizontal  neu¬ 
tral  wind  vector,  for  example,  have  been  binned  and 
averaged  according  to  various  geophysical  indices  in 
order  to  investigate  the  mean  response  of  the  ther¬ 
mospheric  neutral  wind  to  forcings  associated  with 
the  aurora.  A  primary  conclusion  is  that  above  60°  of 
latitude  the  neutral  wind  pattern  in  the  upper  ther¬ 
mosphere  is  best  ordered  in  a  geomagnetic  coordinate 
system  (Killeen  et  al.,  1982 ;  Roble  et  al.,  1983 ;  Hays 
et  al.,  1984;  Rees  et  al.,  198S).  The  high-latitude 
neutral  winds  are  observed  generally  to  follow,  but 
lag  behind,  the  pattern  of  ionospheric  convection, 
with  some  important  differences.  Thus,  in  Fig.  1,  for 
example,  the  neutrals  and  ions  have  qualitatively  simi¬ 
lar  velocities  (speed  and  direction)  throughout  the 
high-latitude  regions,  with  some  important  differ¬ 
ences.  In  general,  for  example,  the  neutrals  have  a 
stronger  sunward  component  in  the  dusk  sector  of 


the  auroral  oval  than  in  the  dawn  sector  due  to  the 
differing  effects  of  the  Coriolis  force  for  anticyclonic 
and  cyclonic  flow  (Killeen  and  Roble,  1984). 

McCormac  et  al.  (1987)  used  several  hundred 
orbits  of  DE  2  over  a  six-month  period  to  investigate 
quantitatively  the  dependence  of  the  high-latitude 
thermospheric  circulation  on  geomagnetic  activity. 
Figure  4  shows  the  measured  average  neutral  wind 
pattern  using  all  available  polar  passes  of  DE  2  for 
active  (A,  >  4)  and  quiet  (K,<  3)  geomagnetic  con¬ 
ditions  in  both  hemispheres.  The  data  are  plotted  in 
geomagnetic  polar  coordinates  (magnetic  local  time 
and  magnetic  latitude).  In  all  four  cases,  the  mean 
neutral  circulation  shows  the  imprint  of  momentum 
transferred  from  the  twin-cell  ionospheric  convection 
pattern,  with  strong  anti-sunward  winds  over  the  geo¬ 
magnetic  polar  cap  bounded  by  strong  sunward  winds 
in  the  dusk  auroral  sector  and  much  weaker  sunward 
or  anti-sunward  winds  in  the  dawn  auroral  sector. 
The  magnitudes  of  the  neutral  winds  are  generally 
greater  for  the  high  case  than  for  the  low  Kp  case, 
as  would  be  expected  from  the  stronger  auroral 
forcings. 

Thayer  et  al.  (1987)  further  investigated  the  sen¬ 
sitivity  of  the  high-latitude  neutral  circulation  pattern 
to  the  sign  of  the  east-west  {B,)  component  of  the 
interplanetary  magnetic  field  (IMF).  They  used  aver¬ 
age  neutral  wind  measurements  from  the  same  data 
base  as  the  McCormac  et  al.  (1987)  study,  discussed 
above,  to  illustrate  changes  in  the  configuration  of  the 
neutral  wind  pattern  that  were  clearly  ordered  by  the 
sign  of  By.  Simultaneous  measurements  of  B,  from  the 
ISEE-3  spacecraft  were  used  to  select  from  and  order 
the  DE  2  neutral  wind  data.  For  this  study,  only  data 
for  which  the  north-south  (F.)  component  of  the  IMF 
was  <  -I- 1  nT  were  used  to  avoid  complications  aris¬ 
ing  from  the  multi-cellular  ion  and  neutral  patterns 
known  to  exist  for  strongly  northward  IMF.  Results 
from  this  study  are  shown  in  Fig.  S,  in  which  data 
from  ~300  orbits  of  DE  2  orbital  passes  have  been 
separated  according  to  the  sign  of  By  and  averaged 
into  bins  of  geomagnetic  latitude  and  geomagnetic 
local  time  for  the  two  hemispheres. 

The  criteria  used  to  separate  the  data  sets  according 
to  By  involved  the  pre-existence  of  a  definite  positive 
or  negative  value  for  By  as  measured  by  ISEE-3  for  1  h 
prior  to  the  DE  orbital  pass.  Passes  occurring  either 
during  very  high  or  very  low  levels  of  geomagnetic 
activity  were  excluded  from  the  study  and  allowance 
was  made  for  the  propagation  time  to  the  Earth  of 
the  IMF  measured  at  ISEE-3  altitudes.  As  can  be  seen 
from  Fig.  S,  there  are  significant  differences  in  the 
mean  thermospheric  circulation  patterns  for  By  posi¬ 
tive  and  negative.  To  some  extent  the  effects  are 


21 


K06 


T.  L.  KiLUiEN  el  al. 


AVERAGED  DE-2  NEUTRAL  WIND  DATA 


(Kp53)  (3+SKp<6) 


Fig.  4.  Averaged  upper  (hermospheric  neutral  wind  measurements  at  ~  .^00 km  altitude  for  (a)  north  pole 
low  K,:  (b)  south  pule  low  At,;  (c)  north  pole  high  <d)  south  pole  high  K,,.  Data  collected  between 
November  IdSI  anu  January  1982  and  between  Novemter  1982  and  January  1983  were  averaged  according 
to  the  given  range  of  and  plotted  in  geomagnetic  polar  coordinates  (magnetic  latitude  and  magnetic 
local  time).  The  outer  circle  of  each  polar  dial  is  at  40  geomagnetic  latitude. 


mirrored  between  hemispheres;  that  is,  a  B,.  positive 
(negative)  signature  in  (he  northern  hemisphere 
resembles  a  B,.  negative  (positive)  signature  in  the 
southern  hemisphere.  In  all  cases,  the  average  cir¬ 
culation  patterns  shown  illustrate  the  dominance  of 
the  dusk  anti-cyclonic  neutral  vortex  mentioned 
above.  The  spatial  magnitude  of  this  vortex,  however, 
is  greater  for  B,.  positive  conditions  in  the  northern 
hemisphere  and  for  B,  negative  conditions  in  the 
southern  hemisphere.  A  distinct  rotation  in  the 
direction  of  the  polar  cap  anti-sunward  Dow  that  was 
dependent  on  the  sign  of  B,.  was  also  noted. 

To  highlight  the  structural  changes  in  the  neutral 
circulation  pattern  that  are  dependent  on  the  sign  of 
the  By  component  of  the  IMF,  Fig.  6  shows  the  aver¬ 
age  anti-sunward  neutral  wind  (i.e.  that  component 
of  the  neutral  wind  parallel  to  the  noon-midnight 
magnetic  meridian)  on  the  nightsidc  of  the  geo¬ 
magnetic  polar  cap  as  a  function  of  magnetic  local 


time  (MLT)  for  the  magnetic  latitude  (MLAT)  range 
72.5-87.5  .  where  the  data  have  been  separated 
according  to  the  sign  of  B,.  for  both  hemispheres 
(Thayer.  1990).  These  plots  illustrate  more  clearly 
the  shift  in  location  of  the  region  of  most  rapid  anti- 
sunward  flow  with  a  change  in  sign  of  B,.  and  dem¬ 
onstrate  the  ‘mirroring'  between  hemispheres  of  this 
characteristic  feature  of  the  neutral  circulation.  In  the 
northern  hemisphere  when  B,.  is  positive  (Fig.  6a),  for 
example,  the  anti-sunward  wind  component  reaches  a 
maximum  in  the  early  morning  sector  at 0400  M  LT. 
When  By  is  negative,  the  maximum  in  the  anti-sun- 
ward  wind  component  occurs  in  the  pre-midnight  sec¬ 
tor  at  ~  2200  MLT.  These  general  results  are  reversed 
in  the  southern  hemisphere.  Fig.  6b.  The  locations  of 
the  various  maxima  in  the  anti-sunward  polar-cap 
neutral  wind  arc  consistent  with  the  locations  of 
maxima  in  the  0, -dependent  convection  electric 
fields  discussed  by  Hkppner  (1977). 
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B  •  DEPENDENCE  OF  DE-2  NEUTRAL  WINDS 
^  (in  geomagnetic  coordinates) 

KpS3 
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Fig.  5.  Averaged  upper  thermospheric  wind  measurements  for  (a)  north  pole  B,  positive :  (b)  north  pole 
B,  negative:  (c)  south  pole  B,  positive:  (d)  south  pole  B,  negative.  Data  collected  during  the  same  period 
as  for  Fig.  4  were  averaged  according  to  the  sign  of  the  B,.  component  of  the  I M  F  and  plotted  in  geomagnetic 
polar  Coordinates  (magnetic  latitude  and  local  time).  The  outer  circle  oP each  polar  dial  is  at  40  geomagnetic 

latitude. 


The  form  of  the  neutral  circulation  for  positive 
values  of  B.  at  high  latitudes  has  also  been  recently 
investigated.  The  existence  of  a  region  of  sunward 
neutral  wind  flow  in  the  central  geomagnetic  polar 
cap  was  noted  by  KtLLEEN  et  at.  (I98S)  for  conditions 
of  persistent  large  positive  values  of  B..  Two  examples 
taken  from  this  study  are  shown  in  Fig.  7.  These  two 
orbits  (7210  and  7211)  occurred  during  a  prolonged 
period  of  strongly  northward  IMF  (B.  ~40nT)  and 
both  the  ion  drift  measurements  (red  bars)  and  neutral 
winds  (blue  arrows)  showed  a  region  of  sunward  flow 
inside  the  geomagnetic  polar  cap.  In  a  later  study, 
McCokmxc'  et  al.  (1991)  systematically  investigated 
the  conditions  for  the  occurrence  of  sunward  polar 
cap  neutral  flow  for  northward  B^  as  a  function  of 
B,  using  the  entire  DE  2  data  base.  They  found  that 
sunward  winds  generally  only  occurred  when  (i)  B. 
was  persistently  positive  and  (ii)  the  magnitude  of  B, 


was  less  than  the  magnitude  of  B..  Figure  8,  taken 
from  the  McCormac  et  al.  study,  is  a  scatter  plot 
illustrating  the  IMF  conditions  for  the  occurrence  of 
sunward  (crosses)  and  anti-sunward  (circles)  neutral 
polar-cap  flow.  A  total  of  66  orbits  were  included  in 
this  study:  in  all  cases  the  hourly  averaged  IMF  B, 
had  been  greater  than  I  nT  positive  for  the  4  h  prior 
to  the  polar-cap  orbital  pass.  It  can  be  seen  from  this 
figure  that  sunward  winds  inside  the  polar  cap  rarely 
occur  when  the  magnitude  of  the  B,  component  is 
greater  than  the  magnitude  of  the  corresponding  B, 
component  (i.c.  most  occurrences  of  sunward  flow 
(crosses)  occur  in  the  uppper  triangular  section  of  the 
plot).  This  complex  dependence  of  the  neutral  wind 
morphology  on  the  relative  magnitudes  of  individual 
IMF  components  is  a  phenomenon  that  clearly  dem¬ 
onstrates  the  strong  coupling  among  the  mag¬ 
netosphere.  ionosphere  and  thermosphere. 
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ANTISOLAR  WIND  COMPONENT  PARALLEL 
TO  NOON-MIDNIGHT  MERIDIAN 
(m/sec) 


Fig.  6,  Average  ami-sunward  neutral  wind  component  par¬ 
allel  to  the  mmn  midnight  meridian  in  meters  per  second  on 
the  night  .side  of  the  geomagnetic  polar  cap  for  both  .signs  of 
B,  in  the  (a)  northern  and  (b)  southern  hemispheres.  The 
wind  components  have  been  averaged  over  magnetic  lati¬ 
tudes  between  72.5  and  117.5  and  for  each  hour  of  magnetic 
local  time  from  IKOO  to  600.  The  figure  is  taken  from  Thavi^r 
(IVW). 


Fig.  K.  Scatter  plot  of  sunward  and  anti-sunward  winds  near 
the  center  of  the  geomagnetic  polar  cap  as  a  function  of  the 
relative  magnitudes  of  the  B,  and  B,.  components  of  the 
interplanetary  magnetic  licld.  A  total  of  66  orbits  are 
included  in  the  plot.  For  all  of  the  included  orbits,  the  hourly 
averaged  IMF  B  had  been  greater  than  I  nT  positive  for  the 
previous  4  h.  This  figure  is  taken  from  MrCoRMAf  el  al. 

(1991). 


2.2.  yorlicity  and  divergence  in  the  high-latitude  ther¬ 
mospheric  wind  pattern 

The  averaged  measurements  of  thermospheric 
winds  from  DE  2  have  enabled  a  study  of  the  vorticity 
and  divergence  in  the  global-scale  mean  neutral  flow 
in  (he  high-latitude  thermosphere  to  be  conducted  for 
(he  first  time  (Thayi-r  and  Killiu:n.  1991).  In  this 
work,  data  from  two  separate  three-month  intervals, 
centered  on  the  winter  solstices  of  1981  1982  and 
1982-198.1.  were  averaged  and  binned  according  to 
geomagnetic  latitude  (MLAT)  and  local  time  (MLT). 
The  DE2  data  were  then  analysed  spectrally  and 
merged  with  VSH  model  predietions  as  described  by 
Thayer  and  Killeen  to  provide  maps  of  the  vorticity 
and  divergence  fields  for  (i)  low  geomagnetic  activity 
and  (ii)  moderately  high  geomagnetic  activity.  These 
maps  are  shown  in  Fig.  9  for  the  northern  hemisphere, 
winter  ease.  Figure  9a  and  c  shows  the  vorticity  maps 
in  geomagnetic  polar  coordinates  for  ^  3  con¬ 
ditions  and  3  4-  ^  $  6  conditions,  respectively, 

and  Fig.  9b  and  d  shows  the  divergence  maps  for  the 
corresponding  eases.  These  fields  have  been  inten¬ 
tionally  smoothed  by  the  spectral  truncation  scheme 
used  in  manipulating  the  satellite  data  and  they  there¬ 
fore  represent  the  mean  vorticity  and  divergence  fields 
for  the  periods  November  1981  -January  1982.  inclus¬ 
ive,  and  November  1982-January  1983.  inclusive. 
Data  from  all  Universal  Times  (UTs)  were  used  to 
generate  these  maps,  based  on  (he  assumption  that 
the  neutral  wind  field  is  largely  invariant  to  UT  when 
plotted  in  geomagnetic  coordinates. 

The  derived  vorticity  patterns  in  units  of  inverse 
seconds,  shown  in  Fig.  9a  and  c,  describe  the  rotation 
in  the  flow  field,  where  negative  (positive)  vorticity 
denotes  anti-cyclonic  (cyclonic)  rotation.  Since  vor¬ 
ticity  is  generated  by  wind  shear,  the  strong  anti- 
sunward  flow  over  the  magnetic  pole,  combined  with 
sunward  flow  in  the  dawn  and  dusk  sectors,  creates 
two  regions  of  counter-rotating  flow.  Negative 
vorticity  or  anti-cyclonic  rotation  occurs  primarily 
in  the  dusk  sector,  with  a  minimum  value  of 

—  518x10  "s  '  for  the  low  K,  case  and 

—  596x10  ‘’s  '  for  the  high  case.  Positive  vorticity 
or  cyclonic  flow  occurs  primarily  in  the  dawn  sector 
with  a  maximum  value  of  426  x  10  ‘s  'for  the  low 
Kp  case  and  559  x  10  "s  '  for  the  high  K,  case.  By 
comparing  (he  vorticity  fields  for  low  and  high  K,.  it 
is  evident  that  an  enhancement  in  geomagnetic  activity 
leads  to  enhanced  vorticity  in  both  dawn  and  dusk 
sectors  and  a  general  growth  in  the  size  of  the  region 
dominated  by  the  twin  vortices.  The  demarcation  line 
of  zero  vorticity  marks  the  region  of  shear  reversals 
where  the  wind  either  reaches  a  maximum  in  the  anti- 
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Fig.  9.  Derived  (a)  vorticity  pattern  for  the  winter  polar  thermosphere  at  ~  300  km  altitude  for  conditions 
of  low  geomagnetic  activity  (northern  hemisphere) ;  (b)  corresponding  divergence  pattern ;  (c)  vorticity 
pattern  for  moderately  active  geomagnetic  conditions ;  (d)  corresponding  divergence  pattern.  The  patterns 
are  shown  in  geomagnetic  polar  coordinates  (magnetic  latitude ;  pole  to  40"  and  magnetic  local  time).  The 
contour  interval  is  lOOx  l()"''s''  for  the  vorticity  plots  and  50x  I0“‘s"'  for  the  divergence  plots.  These 
plots  were  generated  from  spectrally  analysed  OE2  data  and  VSH  model  predictions  (from  Thavkr  and 

Killeen.  1991). 


sunward  direction  in  the  central  polar  cap.  or  a 
maximum  in  the  sunward  direction  in  the  dawn  and 
dusk  sectors.  Thus,  the  transitions  in  the  vorticity 
fields  may  be  used  to  define  the  spatial  dimensions  of 
the  high-latitude  neutral  circulation  pattern.  Below 
~60  ML  AT,  the  vorticity  values  are  relatively  small. 
The  dominance  of  the  signatures  of  the  twin  vortex 
motion  induced  in  the  high-latitude  thermosphere  by 
ion  drag  processes  is  evident  in  these  maps  which 
are  strongly  reminiscent  of  similar-looking  maps  of 
equipotentials  for  the  polar  electric  fields  (e.g. 
Heppner.  1977;  Heelis  er  a/.,  1982). 


The  derived  divergence  patterns,  shown  in  Fig.  9b 
and  d.  describe  the  rate  of  change  of  area  in  the  flow 
fields  determined  by  the  longitudinal  difference  in  the 
zonal  wind  components  and  the  latitudinal  difference 
in  the  meridonal  wind  components,  where  negative 
values  denote  convergence  in  the  flow  field  and  posi¬ 
tive  values  denote  divergence.  Once  again,  these  fields 
represent  mean  patterns  for  the  epoch  in  question  and 
are  plotted  in  units  of  inverse  seconds.  It  can  be  seen 
that  the  divergence  values  are  generally  smaller  (by 
about  a  factor  of  4)  than  the  vorticity  values  discussed 
above.  This  is  due  to  the  greater  importance  of  the 
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largely  divergence-free  ton-drag  momentum  term  in 
establishing  the  neutral  circulation.  Perturbations  in 
the  divergence  fields  arc  primarily  forced  by  localized 
regions  of  heating,  which  can  lead  to  upwclling  and 
positive  divergence,  or  regions  of  relative  cooling, 
which  can  lead  to  downwelling  or  negative  divergence 
(convergence) .  The  divergence  patterns  are  more  com¬ 
plex  than  their  vorticity  counterparts  and  are  char¬ 
acterized  by  localized  islands  of  positive  and  negative 
divergence  associated  with  regions  of  high-latitude 
heating  and  subsidence.  The  effect  of  enhanced  icvcis 
of  geomagnetic  activity  is  to  widen  the  area  of  large 
perturbations  in  the  divergence  fields  and  to  intensify 
the  maxima  and  minima  observed.  In  these  localized 
regions,  the  divergence  field  becomes  commensurate 
with  the  vorticity  field  in  magnitude. 

The  ability  to  generate  thermospheric  vorticity  and 
divergence  maps  such  as  those  shown  in  Fig.  9.  pro¬ 
vides  a  new  tool  for  the  interpretation  of  high-latitude 
thermospheric  circulation  patterns,  allowing  for  a 
more  detailed  kinematic  analysis  of  the  flow  con¬ 
figurations.  In  particular,  this  novel  tool  gives 
additional  insight  into  the  relative  strengths  of  the 
various  sources  of  momentum  and  energy  responsible 
for  driving  the  winds. 


3.  TMF.  VECTOR  SPHERICAL  HARMONIC  MODEL 

The  Vector  Spherical  Harmonic  (VSH)  model  has 
been  described  by  Killef.n  et  al.  (1987).  It  is  a  semi- 
empirical  'hybrid'  model  that  is  based  on  numerical 
simulations  of  the  NCAR-TICCM  and  available 
experimental  data  from  DE2.  The  NCAR-TIGCM 
solves  the  primitive  equations  of  dynamic  meteor¬ 
ology  adapted  to  thermospheric  altitudes,  including 
equations  for  the  conservation  of  energy  and  momen¬ 
tum  and  the  individual  species  continuity  equations. 
The  physics  inherent  in  the  model's  parameterizations 
and  input  prescriptions  arc  those  appropriate  to  tlwr- 
mospheric  altitudes.  The  basic  model  and  subsequent 
major  developments  have  been  described  in  detail  by 
Dickinson  et  al.  (1981,  1984),  Roble  ei  al.  (1982). 
Fesen  el  al.  (1986)  and  Roble  and  Ridley  (1987). 
Most  recently,  the  model  was  further  extended  (Roble 
et  al.,  1988)  to  include  a  self-consistent  aeronomic 
scheme  for  the  thermosphere  and  ionosphere  using  a 
fully  Eulerian  approach. 

The  VSH  model  uses  output  from  various  'generic' 
runs  of  the  TIGCM  to  develop  a  library  of  truncated 
spectral  coefficients  that  represent  the  global  ther¬ 
mospheric  wind,  temperature  and  neutral  com¬ 
position  fields  for  a  set  of  diflerent  geophysical  con¬ 
ditions  (solar  maximum,  solar  minimum,  equinox. 


solstice,  etc.).  The  three-dimensional  wind  field  is  rep¬ 
resented  by  a  vector  spherical  harmonic  expansion  in 
the  horizontal,  a  Fourier  expansion  in  Universal  Time 
and  a  polynomial  expansion  in  altitude.  The  global 
temperature  field  differs  in  that  a  scalar  spherical  har¬ 
monic  expansion  is  used  in  the  horizontal  and  a  Bales 
model  profile  is  used  in  the  vertieal.  Neutral  com¬ 
position  fields  are  described  using  logarithmic  scalar 
harmonic  expansions.  The  thermospheric  state  vari¬ 
ables  (wind,  temperature  and  composition)  may  then 
be  calculated  simply  using  a  computer  subroutine 
which  reads  in  the  coefficients  and  regenerates  the 
required  field  as  a  function  of  space  and  time  and 
geophysical  condition  (^10.7  value  and  A^).  Recently, 
data  from  DE  2  have  also  been  incorporated  into  the 
VSH  model  description  using  an  objective  analysis 
scheme  to  merge  the  experimental  data  and  theor¬ 
etically  calculated  TIGCM  gridded  output  fields,  in 
this  objective  analysis  scheme,  measurements  from 
DE  2  are  used  whenever  available,  and  TIGCM  model 
results  are  used  in  the  absence  of  experimental  data, 
with  the  whole  process  being  carried  out  in  the  spectral 
domain  to  provide  for  smooth  merger  of  the  exper¬ 
imental  and  theoretical  results. 

3. 1 .  VSH  model  comparisons  with  thermospheric  wind 
observations 

An  example  of  a  comparison  between  the  global 
output  of  the  VSH  model  and  averaged  DE  2  data  is 
shown  in  Fig.  10.  The  lower  part  of  this  figure  (lOb) 
shows  the  bin-averaged  neutral  wind  field  from  FPI 
and  WATS  as  a  function  of  UT  and  latitude,  coded 
as  vectors,  for  the  same  two  three-month  study 
periods  as  for  the  previous  results  shown  earlier.  The 
thermospheric  temperatures  from  FPI  and  WATS  are 
coded  as  colors  according  to  the  scale  at  the  right.  As 
can  be  seen,  the  largest  wind  speeds  are  observed  in 
the  two  high-latitude  regions  (southern  hemisphere 
polar  regions  towards  the  lower  part  of  the  plot  and 
northern  hemisphere  polar  regions  towards  the  upper 
part).  The  measured  thermal  structure  indicates  a  rela¬ 
tively  hot  summer  polar  region  (temperatures  up  to 
~  I SOO  K)  and  a  relatively  cold  low-latitude  region 
(blue  band  in  middle  of  plot — temperatures  ~900  K). 
For  comparison.  Fig.  I  Oa  shows  the  equivalent  output 
from  the  VSH  model,  where  the  model  was  interro¬ 
gated  to  provide  calculated  December  solstice  wind  and 
temperature  fields  in  the  same  format  as  the  exper¬ 
imental  data.  As  can  be  seen,  there  is  reasonable  agree¬ 
ment  between  the  model  and  the  global-scale  exper¬ 
imental  data.  The  seasonal  thermal  structure  (sum- 
mer-to-wintcr  gradient)  is  present  in  the  model  which, 
however,  tends  to  underestimate  the  northern  polar 
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temperatures.  The  wind  pattern  shows  reasonable 
agreement  also,  though  there  is  considerably  more 
structure  in  the  experimental  results  than  in  the  model 
predictions.  The  morphology  or  the  high>latitudc 
flows  and  reversals  arc  generally  well  matched  by  the 
model. 

To  provide  additional  comparisons  with  the  new 
VSH  model.  Figs  1 1  and  12  show  comparisons  with 
ground-based  thermospheric  neutral  wind  measure¬ 
ments.  Fig.  I  la  shows  thermospheric  wind  measure¬ 
ments  from  two  optical  observatories  in  Greenland, 
one  at  Thule  (76.0  N,  70.0  W)  in  the  central  region 
of  the  geomagnetic  polar  cap  and  one  at  Sondrestrom 
(67.0  N,  51.0  W)  on  the  northern  edge  of  the  auroral 
aonc.  The  winds  shown  represent  bin-averaged  data 
for  geomagnetic  quiet  times  {K^  <  .1),  taken  during 
the  1987  1988  observing  season  (September- March). 
The  vectors  shown  represent  bin-averaged  winds  at 
^250  km  altitude  from  both  stations  for  all  available 
clear  nights  of  observations  during  low  geomagnetic 
activity  conditions  {K^  <  3),  plotted  in  polar  geo¬ 
graphic  coordinates.  The  inner  set  of  vectors  show  the 
averaged  winds  from  Thule  and  the  outer  set  of  vec¬ 
tors  show  the  winds  from  Sondrestrom.  For  compari- 
.son.  Figure  1 1  b  shows  the  output  from  the  VSH  model 
which  was  interrogated  to  match  the  conditions 
appropriate  to  the  experimental  results.  The  measured 
winds  show  a  largely  anti-sunward  flow  over  the  polar 
cap,  bounded  on  the  dusk  side  with  a  region  of  sun¬ 
ward  flow,  and  with  a  turning  and  abatement  of  the 
winds  on  the  dawn  side  of  the  polar  cap.  The  com¬ 
parison  between  experiment  and  model  indicates 
reasonable  agreement,  with  the  major  characteristics 
of  the  flow  pattern  being  well  represented  by  the 
model.  There  are  some  discrepancies,  however, 
notably  in  the  direction  of  the  strong  anti-sunward 
flow  over  Sondrestrom  near  midnight. 

Figure  12  shows  experimental  neutral  wind  data 
from  the  incoherent  scatter  radar  at  Sondrestrom, 
separated  aocordtAg  to  the  sign  of  the  IMF  0,.  com¬ 
ponent.  These  data  represent  approximately  2200  h  of 
measurements  obtained  between  April  1983  and  July 
1988.  The  radar  winds  are  averaged  between  210  and 
360  km  altitude.  The  VSH  model  was  interrogated  to 
provide  predictions  for  the  neutral  winds  during  these 
periods  and  these  results  are  shown  in  Fig.  12.  The 
agreement  between  the  VSH  model  predictions  and 
the  radar-determined  winds  for  the  given  site  is  again 
reasonable,  perhaps  with  the  exception  of  the  early 
morning  period  for  0,  positive  where  the  model  over¬ 
estimates  the  magnitude  of  the  measured  wind.  Since 
the  empirical  component  of  the  VSH  model  is  pro¬ 
vided  by  DE2  measurements  referred  to  the  1981- 
1983  solar  maximum  period,  it  is  encouraging  to  note 
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Pig.  II.  (a)  Measured  neutral  winds  at  ~ 250 km  altitude 
determined  from  the  Fabry- Perot  observatories  at  Thule, 
and  Sondrestrom.  Greenland,  during  the  observing  season 
of  1987  1988.  The  data  represent  bin-averaged  winds  for  all 
ctear-sky  observations  from  the  two  stations.  The  winds  are 
coded  as  arrows  (scale  at  lower  right)  with  the  base  of  the 
arrow  located  at  the  appropriate  local  time  and  latitude  for 
the  particular  station.  The  Thule  measurements  constitute 
the  inner  circle  of  vectors  and  the  Sondrestrom  measure¬ 
ments  constitute  the  outer  circle  of  vectors,  (b)  VSH  model 
predictions  for  the  corresponding  geophysical  conditions. 


the  degree  of  model  agreement  with  neutral  wind  data 
obtained  from  a  quite  different  time  period  and  phase 
of  the  solar  cycle. 

4.  THERMOSPHKRIC  COMPOSITION  PERTURBATIONS 
AT  HICH  LATITUDES 

The  resulls  described  above  have  demonstrated 
various  characteristics  of  the  mean  thermospheric 
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Fig.  12.  Comparison  of  Sondrestrom  incoherent  scatter  radar  average  meridional  winds  (in  the  magnetic 
meridian)  with  predictions  of  the  Vector  Spherical  Harmonic  Model  (Kiluxn  et  a!.,  1987)  as  a  function 
of  IMF  0,  orientation  for  B.  negative  conditions.  The  radar  data  set  includes  approximately  22(X)h  of 
measurements  obtained  from  April  1983  to  July  1988.  Radar  winds  were  (irst  averagni  over  hourly  intervals 
from  210  to  360  km  for  each  experiment.  Theaverage  winds  as  a  function  of  UT  and  IMF  orientation  were 
then  determined  by  binning  the  entire  data  set  as  a  function  of  UT  hour  and  average  IMF  orientation 
during  the  period  from  I  to  2  h  prior  to  the  beginning  of  each  hourly  bin.  The  figure  is  taken  from  Johnson 

ere/.  (1991). 
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neutral  circulation  at  high  geomagnetic  latitudes.  The 
wind  patterns  show  a  circulation  that  has  a  strong 
rotational  component  associated  with  the  twin-cell 
ionospheric  convection  pattern,  with  a  significant 
divergent  component  driven  largely  by  pressure  gradi¬ 
ent  forcing,  it  might  be  expected  that  the  large  mag¬ 
nitude  horizontal  neutral  winds  and  associated  ver¬ 
tical  motions  would  have  a  significant  effect  on  the 
compositional  balance  of  the  upper  thermosphere.  To 
investigate  this  possibility.  Burns  et  al.  (1989)  have 
used  the  NCAR-TIGCM  to  study  the  changes  in  ther¬ 
mospheric  composition  that  result  from  wind  systems 
such  as  discussed  above.  Figure  13  shows  an  example 
of  a  set  of  calculations  from  their  work,  illustrating  the 
various  dynamical  forcing  terms  for  thermospheric 
composition  variations  during  and  in  the  aftermath 
of  a  simulated  geomagnetic  storm. 

This  figure  shows  NCAR-TIGCM  calculations  of 
the  terms  responsible  for  changes  in  the  mass  mixing 
ratio  of  N;  (Vm,)  during  the  large  geomagnetic  storm 
that  occurred  on  24  November  1982.  The  terms  are 
referred  to  an  altitude  of  ^  280  km  and  the  location 
of  the  geomagnetic  pole  is  shown  as  the  cross  in  Fig. 
13d.  The  plots  show  the  differences  between  storm 
and  quiet  times  for  the  following  terms :  (a)  vertical 
advcctive  forcing  of  (b)  horizontal  advective 
forcing  of  yN. :  (c)  molecular  diffusion  effects  and  (d) 
yN.  itself.  The  first  three  plots  (I3a-c),  showing  com¬ 


positional  adjustment  forcing  terms,  are  for  14  UT, 
some  3  h  after  the  comnnencement  of  the  storm,  and 
the  plot  showing  the  mass  mixing  ratio  (13d)  is  for 
19  UT  or  2  h  into  the  storm  simulation.  Thus,  these 
plots  illustrate  some  of  the  physical  forcing  processes 
that  give  rise  to  the  changes  in  mass  mixing  ratio  at 
~  280  km  a  few  hours  after  the  commencement  of  a 
large  geomagnetic  storm. 

The  pattern  of  compositional  forcing  that  is  seen  at 
high  latitudes  for  this  storm-time  case  results  from 
both  horizontal  and  vertical  advection,  where  the 
horizontal  and  vertical  winds  are  due  primarily  to 
ion  drag  and  Joule  heating,  respectively.  Figure  13a 
illustrates  the  rate  of  change  of  mass  mixing  ratio  due 
to  the  advective  effects  of  vertical  winds.  A  region  of 
positive  values  (associated  with  upwelling  near 
regions  of  enhanced  heating)  can  be  seen  to  be  roughly 
coincident  with  the  location  of  the  aurora)  zone  and 
regions  of  negative  forcing  (associated  with  sub¬ 
sidence)  can  be  seen  to  lie  equatorward  of  the  auroral 
zone  and  within  the  geomagnetic  polar  cap.  Figure  13b 
illustrates  the  rate  of  change  of  y^,  due  to  the  advec¬ 
tive  effects  of  horizontal  winds.  The  region  of  large 
positive  values  in  the  post-midnight  sub-auroral  sec¬ 
tor  is  due  to  the  existence  of  the  post-midnight  surge 
in  the  meridional  winds  which  tends  to  increase  the 
mass  mixing  ratio  of  nitrogen  downstream  from  the 
polar  regions.  Figure  1 3c  illustrates  the  rate  of  change 
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Fig.  t  NCAR-TICiCM  vakulalions  oF  ihc  forcing  terms  responsible  for  changes  in  the  mass  mixing  ratio 
of  N<  (*I'k.)  during  Ihc  large  geomagnetic  storm  that  occurred  on  24  November  19X2.  The  plots  show  the 
mas.vmixing-ratio  forcing-term  dilTcrenccs  between  storm  and  quiet  times  for  the  following :  (a)  vertical 
advcction;  (b)  horizontal  advcclion:  (c)  molecular  difru.sion  and  (d)  (H'n.).  The  lirst  three  terms  are  for 
14  UT.  some  .1  h  after  the  commencement  of  Ihc  storm,  and  Ihc  last  plot,  showing  mass  mixing  ratio,  is  for 
19  UT.  The  cross  in  Fig.  I  .kl  indicates  the  kwalion  of  the  geomagnetic  pole. 


of  due  to  molecular  dilTusion.  This  term  lends  to 
be  negative  throughout  (he  high-la(i(ude  region  as  i( 
attempts  to  restore  (he  thermosphere  to  the  ditfusivc 
equilibrium  that  has  been  interrupted  by  the  high- 
latitude  dynamical  processes.  Maximum  values  of  (he 
perturbed  H'n.  (Fig.  I.^d)  are  found  in  (he  auroral  and 
sub-auroral  morning  sector,  in  agreement  with  (he 
DE  2  observations  reported  by  Hkdin  and  Carkinan 
(l9g.S).  The  magnitudes  of  the  three  terms  shown  in 


Fig.  I.^a-c  arc  roughly  commensurate,  indicating  that 
all  three  proces.ses  (vertical  and  horizontal  advection 
and  molecular  diffusion)  arc  of  roughly  equal  import¬ 
ance  in  (he  establishment  of  the  storm-time  com¬ 
positional  state  of  (he  high-latitude  thermosphere. 
This  type  of  quantitative  compositional  forcing  term 
analysis  thus  enables  a  deeper  insight  into  the  coupling 
of  dynamics  and  composition  in  (he  high-latitude 
thermosphere. 
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5.  CONCLUDING  RrJMABKS 

Recent  experimental  results  from  the  Dynamics 
Explorer-2  spacecraft  and  ground-based  optical  and 
radar  observatories  have  illustrated  the  manner  by 
which  the  high-latitude  processes  of  ion  drag  and 
Joule  heating  control  the  behavior  of  the  high-latitude 
upper  thermosphere.  The  neutral  winds  at  high  alti¬ 
tudes  typically  lend  to  follow  the  twin-cell  pattern  of 
ionospheric  convection,  as  momentum  is  transferred 
from  the  ions  to  the  neutrals  at  a  rate  that  is  pro¬ 
portional  to  the  ion  density.  As  a  consequence  of  this 
tight  ion-neutral  momentum  coupling,  the  neutral 
wind  pattern  displays  strong  dependencies  on  the  level 
of  geomagnetic  activity  and  on  the  orientation  of  the 
interplanetary  magnetic  held,  much  in  the  same  way  as 
the  ionospheric  convection  pattern  does  itself.  Neutral 
wind  vortices  are  established  whose  geometries  and 
intensities  depend  on  the  nature  and  time-history  of 
the  momentum  forcing.  Upper  thermospheric  tem¬ 
peratures  are  directly  controlled  by  the  Joule  and  par¬ 
ticle  heat  sources  associated  with  auroral  latitudes. 
The  auroral  heating  also  serves  to  modulate  the  ther¬ 
mospheric  circulation  through  the  generation  of 


strong  pressure-gradient  forces  and  through  local 
upwelling  of  the  atmosphere.  Changes  in  the  vertical 
and  horizontal  winds  at  high  latitudes  lead,  in  turn, 
to  neutral  compositional  perturbations.  The  various 
dynamical,  energetic  and  compositional  variations 
associated  with  auroral  processes  are  thus  highly 
coupled  with  each  other. 

Significant  progress  has  been  made  recently  in  the 
development  of  semi-empirical  models  of  the  ther¬ 
mosphere  that  utilize  the  large  and  comprehensive 
new  data  resources.  One  of  these  models,  the  Vector 
Spherical  Harmonic  model,  has  been  shown  to  pro¬ 
vide  reasonably  accurate  predictions  for  global-scale 
thermospheric  wind  systems. 
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ABSTRACT 

Ground  based  incoheieat  scatter  radar  (ISR)  and  Fidiry  Perot  interferometer  (FPI)  shidies  in  the  northern 
high  ladtades  during  the  period  1983  to  1989  have  ^own  that  the  F-region  neutral  wind  field  pattern 
dqpends  upon  the  sign  of  the  IMF  parameters.  Fbr  example,  the  cell  structure  of  the  northern  hemisphere 
hi^  latitude  neutral  wind  field  during  periods  of  low  geomagnedc  activity  depends  to  a  large  degree  iq>on 
the  sign  of  the  IMF  By  parameter. 

Long  term  monitoring  of  the  F-tegion  thermosphere  by  FPI  in  Thule,  Greenland,  and  by  both  FPI  and 
ISR  in  Spndre  StrpmQord,  Greenli^,  have  made  it  possible  to  produce  nups  of  average  meridional  and 
zonal  wind  fields  for  various  IMF  configurations  for  northern  high  latitudes.  Comparison  of  observations 
with  theoretical  wind  field  modelling,  such  as  the  Vector  Spherical  Harmonic  mo^l,  indicates  that  most 
observed  features  are  consistent  with  the  models. 

INTRODUCTION 

Neutral  winds  in  the  polar  ctq)  F-region  are  controlled  by  a  combination  of  pressure  gradient  forces  and 
ion  drag  fortes,  the  laber  being  the  primary  driving  force  at  high  latitudes.  Coupling  biRween  the  ion  and 
the  neutral  constituents  in  the  F-tegion  thetmospheie,  which  is  dependant  upon  the  ion  neutral  collision 
ftequeiKy,  fortes  the  neutral  component  to  follow  the  ion  drift  convection  pattern  with  an  appropriate 
time  constant  As  a  result  the  F-tegion  neutral  wind  circulation  pattern  for  high  latitudes  usually  exhibits 
a  similarity  to  the  ion  convection  cell  pattern. 

The  configuntio/i  of  the  Interplanetary  Magnetic  Reid  (IMF)  has  a  iiuuked  effect  upon  the  size  and 
geometry  of  the  ion  convection  ceil  pattern,  and  consequently,  the  neutral  circulation  pattern.  In 
particular,  the  sign  of  the  By  component  of  die  IMF  controls  the  relative  size  of  the  dawn  and  dusk 
ionospheric  circulation  cells.  However,  the  time  constant  for  momentum  transfer  Bom  ions  to  neutrals  is 
of  the  order  of  several  hours,  implying  that  rapid  changes  in  the  IMF  are  not  immediately  apparent  in  the 
neutral  wind  pattern.  Several  studies  have  been  published  which  document  the  influence  of  the  By 
component  on  the  high  latitude  thermoqihere:  ground-based  optical  /1-S/,  ground-based  radar  /b/,  satellite 
m,  and  theoretical  srodies  /8/. 

In  this  study,  data  acquired  by  tqitical  techniques  at  S^idre  Strpmfiord  and  Thule,  Greenland  are  used  to 
construct  averaged  horizontal  winds  for  both  signs  of  the  IMF  By  component  These  results  are  compared 
with  predicted  thermospheric  circulation  patterns  from  a  Vector  Spheral  Harmonic  (VSH)  model.  The 
comparison  indicates  that  this  model  provides  a  reasonable  representation  of  the  optically  observed 
horizontal  wind  components  in  the  thennosphere  for  the  two  signs  of  the  IMF  By  component  Fitully,  the 
averaged  neutral  winds  derived  from  the  incoherent  scatter  radar  at  Spndre  Strpntfjotd  for  the  same 
geomagnetic  conditions  are  compared  with  the  optical  measurements  and  the  modelled  winds. 

RESULTS 


Similar  Fabry  Perot  interferometers  are  located  in  both  Sendre  StrpmQotd  and  Thule,  Greenland,  the 
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fonner  collocated  with  an  incoherent  scatter  radar  system.  The  optical  instrument  is  described  in 
Meriwether  et  al.  14/.  The  interferometers  were  in  routine,  automated  operation  during  the  solar 
minimum  period  acquiring  F-region  neutral  wind  measurements  during  dark  sky  periods  through  the 
entiK  (^Mical  observing  season  -  September  to  April.  Goud  cover  dau  for  both  stations  was  obtained  from 
the  Danish  Meteorological  Institute  and  the  United  States  Air  Fbrce  for  Spndie  StrpmQord  and  Thule, 
reflectively.  This  information  permitted  the  removal  of  data  acquired  during  non-ideal  observing  periods. 
Geomagnetic  indices  and  IMF  values  (in  Geocentric  Solar  Magnetofiheric  coordinates)  were  obtained 
from  the  NSSDC  data  base.  In  total,  four  years  of  Spndre  Strpmfjord  data  and  three  years  of  Thule  data 
during  solar  minimum  were  used  to  generate  the  finaJ  averages.  Ilie  horizontal  wind  components  may  be 
considered  as  representative  of  the  neutral  wind  pattern  for  the  northern  polar  cap  for  conditions  of 
average  geomagnetic  activity  during  the  solar  minimum  period. 

The  conditions  chosen  for  the  binning  process  are  similar  to  those  used  by  Thayer  et  al.  P/.  The  IMF 
parameters  used  in  the  current  sorting  schetrte  were  one  hour  averages  for  both  By  and  B,,  which  is  the 
same  resolution  as  the  cloud  cover  data.  Each  hwizontal  wind  measurement  was  assigned  By,  B,,  K^,  and 
cloud  index  values.  The  winds  were  then  sorted  into  one  hour  wide  l»ns  based  upon  the  following  criteria: 
a)  low  geomagnetic  activity  (Kp  ^  3);  b)  southward  Bj  component  (Bj  <  1  nT);  c)  the  sign  of  By  constant 
for  at  least  two  hours  preceding  the  wind  tneasur«nent;  d)  no  inclusion  of  wind  data  within  the  bins  if 
there  ate  no  IMF  data  available  for  the  measurement  period;  e)  cloud  cover  must  be  low.  This  strict 
filtering  criteria  reduced  the  neutral  wind  data  set  to  slightly  over  2%  of  its  original  size. 

figures  la)  and  lb)  display  the  results  of  the  sorting  analysis.  The  figures  are  shown  using  a  geontagnetic 
latitudeAnagnetic  local  time  coordinate  system  in  a  polar  dial  display.  The  length  of  the  vector  is 
proportional  to  the  magnitude  of  the  horizontal  wind  as  indicated  by  the  vector  in  tte  bottom  right-hand 
comer.  The  tail  of  a  vector  is  located  at  the  geonugnetk  latitude  of  the  station  and  the  magnetic  hour  of 
the  bin.  The  outer  ring  of  vectors  c^'ie^ponds  to  F-region  neutral  winds  over  Spndre  StrpmQord,  while 
the  inner  ring  shows  the  same  for  Thule.  Even  though  Spndre  Strpmfjord  is  at  the  arctic  circle,  there  are 
still  periods  at  winter  solstice  during  which  the  sky  illumination  is  too  high  to  perform  data  acquisition.  As 
a  result,  this  ring  of  data  does  not  cover  an  entire  day.  Figure  la)  refen  to  By  negative  conditions,  while 
Rgure  lb)  shows  By  positive  results. 


s«u<lri  Slromriard  ■u4  ThuU  rpl  —  »ynn*|)4«  Sondrt  SUremlJord  anii  Tbult  FPI  “  *ypp*P'» 
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Fig.  1.  a)  Polar  dial  plot  representing  average  horizontal  neutral  wind  vectors  in  the  F- 
region  thermosphere.  The  coordinate  system  is  geomagnetic  latitude/corrected  geomagnetic 
local  time.  The  plot  is  centered  on  the  north  geomagnetic  i»lc  with  circles  shown  for  80® 
and  70°  latitude.  The  magnitude  of  the  horizontal  neutral  wind  is  proportional  to  the  length 
of  the  vector,  as  shown  by  the  arrow  in  the  lower  right  comer.  The  location  of  the  tails  of 
the  vectors  correspond  to  the  geomagnetic  latitudes  of  Thule  (inner  circle)  and  Snndre 
StrpmQord  (outer  circle).  This  panel  is  for  By  negative  conditions,  filtered  as  described  in 
the  text  b)  Same  as  n)  but  for  By  positive  conditions. 

The  neutral  wind  model  chosen  for  the  comparison  was  the  VSH  model  /9/.  This  is  a  computer  model 
which  provides  time  dependent,  horizontal  neutral  wind  vectors.  This  model  consists  of  a  set  of 
coefficients  derived  from  geophysical  field  outputs  generated  by  a  selection  of  various  NCAR 
thermospheric  general  circulation  model  runs.  A  vector  spherical  harmonic  expansion  of  the  coefficients 
is  used  to  represent  the  wind  field,  while  a  Fourier  expansion  represents  the  temporal  dimension. 
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In  this  study,  two  diumatly  reproducible  NCAR  TGCM  “steady  state"  runs  were  chosen  for  comparison 
with  the  optically  measur^  winds,  one  for  the  By  positive  case  and  the  other  for  the  By  negative  case. 
Both  runs  were  made  for  conditions  of  average  geomagnetic  during  the  solstice  periods.  Average 
geomagnetic  activity  was  defined  as:  a)  Kp=  3;  b)  Ap  =  20;  c)  cross-cap  potential  =  30  kV.  In  both  runs, 
the  absolute  magnitude  of  By  was  7  nT. 

Neutral  winds  derived  from  incoherent  scatter  radar  dau  for  Sondre  Stremfjord  have  previously  been 
presented  for  similar  conditions  by  de  la  Beaujardiire  and  Wickwar  Ibl.  In  their  study,  neutral  winds  from  a 
two  year  set  were  used.  In  the  current  study,  neutral  winds  derived  from  the  entire  radar  data  base  (April, 
^3  to  July,  1988)  were  used  to  construct  averaged  neutral  winds  in  the  geomagnetic  meridional  direction’. 
Prior  to  sorting  of  the  winds  according  to  IMF  conditions,  the  winds  for  range  gates  between  210  and  360  km 
were  averaged  to  generate  a  single  value  representing  the  neutral  thermospheric  meridional  wind  in  the  F- 
r^ion.  Next,  the  winds  were  sorted  into  one  hour  wide  bins  using  similar  By  and  B*  criteria  as  for  the  nprical 
data  sotting,  the  only  difference  being  that  the  hourly  averaged  wind  measurements  were  associated  with  the 
average  of  the  IMF  parameters  obtained  1  to  2  hours  prior  to  the  wind  measurements. 

A  comparison  of  the  observed  neutral  winds  with  the  modeled  neutral  winds  is  shown  in  Figures  2a)  and  2b), 
the  former  for  By  negative  conditions  and  the  latter  for  By  positive  conditions.  Here,  only  the  averaged  geo¬ 
magnetic  meridional  winds  as  a  function  of  time  for  Spndre  SirnmQoid  are  shown.  The  srarisri<-a|  uncertainty 
in  the  mean  is  also  indicated  for  both  the  radar  (R)  and  optical  (O)  averages.  There  are  no  optical  measure¬ 
ments  between  13  and  18  hours  Universal  Time  since  this  period  brackets  local  noon  in  S0ndre  Sti0mQord. 


geomagnetic  meridian.  The  average  neutral  winds  acquired  experimentally  are  shown  with 
their  respective  statistical  uncertainties  in  their  means;  O  refers  to  FPJ  data  and  R  refers  to 
ISR  data.  The  curve  displays  the  corresponding  model  winds  from  the  VSH  code.  This 
panel  refers  to  By  negative  conditions,  b)  Same  as  a)  but  for  By  positive  conditions. 

DISCUSSION 

Several  features  are  obvious  from  the  polar  dial  plots  of  Rgure  1  and  the  collated  data  in  Figure  2.  For 
By  positive  conditions,  the  magnitude  of  the  horizontal  neutral  winds  is  smaller  in  Thule  than  for  By 
negative  periods.  In  addition,  the  magnitude  of  the  meridional  winds  measured  at  Sondre  Strpmfjord  is 
smaller  for  By  positive.  Although  not  shown  here,  the  absolute  magnitude  of  the  experimental  averaged 
zonal  winds  for  Spndre  Stromfjotd  is  larger  for  the  By  positive  case,  unlike  the  meridional  averages. 
Assuming  a  two  cell  structure  in  the  neutral  horizontal  wind  flow  within  the  polar  cap  and  fitting  that  with 
the  observed  horizonul  wind  vectors  indicates  that  during  the  By  negative  case  the  cells  appear  to  be 
similar  in  size.  Likewise,  for  the  By  positive  condition,  an  asymmetry  is  apparent  between  the  two  cells, 
the  dusk  cell  being  much  larger  in  size  than  the  dawn  cell,  encroaching  into  the  dawn  sector. 

The  behavior  of  the  presumed  cell  structure  of  the  neutral  wind  field  as  a  function  of  the  sign  of  IMF  By 
within  this  study  has  some  similarities  with  the  comprehensive  study  performed  by  Thayer  et  at.  PI.  Their 
work  is  based  upon  DE2  satellite  measurements  of  the  neutral  wind  field  during  solar  maximum.  The 
asymmetry  in  the  cell  pattern  during  By  positive  in  the  northern  hemisphere  observed  by  the  satellite  is 
also  evident  during  the  current  study's  solar  minimum  observations.  They  also  observed  a  symmetric 
pattern  during  By  negative  conditions.  However,  the  large  surge  that  they  observed  in  the  magnitude  of  the 
winds  at  ~7S°  invariant  latitude  during  By  positive  periods  is  not  evident  during  the  solar  minimum  study. 
Similarities  between  the  averaged  wi<v)  vectors  in  Figure  la)  and  lb)  and  the  data  presented  by 
Meriwether  and  Shih  PI,  for  Spndre  Stromfjord  and  Meriwether  et  al.  /4/,  for  Thule  are  also  evident.  In 
the  latter  two  studies,  individual  nights  of  data  were  shown  to  display  a  By  dependence:  an  enlargement  in 
the  evening  cell  towards  the  morning  cell  for  By  positive  conditions. 
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The  observed  geomagnetic  meridional  neutral  winds  shown  in  Figures  2a)  and  2b)  show  reasonable 
agreement  with  the  VSH  model  winds  during  dark  conditions.  In  both  the  By  negative  and  the  By  positive 
cases,  the  phase  of  the  observed  wind  is  similar  to  the  modelled  wind  with  zero  crossings  differing  by  no 
more  than  about  one  hour.  The  magnitude  of  the  observed  winds  during  dark  hours  is  similar  to  the 
modelled  winds,  but  both  the  optical  and  the  radar  winds  show  structure  in  the  binned  averages.  The  fact  that 
the  optical  winds  differ  from  the  radar  neutral  winds  is  not  surprising  since  there  were  very  few,  if  any, 
coordinated  measurements  between  the  two  sets  of  observations.  The  greatest  discrepancy  between  the 
model  and  the  observations  occurs  during  daylight  conditions.  This  also  is  not  too  surprising  when 
consideration  is  made  regarding  the  geometry  of  the  auroral  oval  with  respect  to  Spndre  Strpmfj^.  The 
cusp/cleft  region  passes  through  the  geomagnetic  meridian  over  this  station  at  approximately  1400  UT 
which  corresponds  to  the  period  of  greatest  discrepancy.  The  modelled  results  provide  wind  vectors  on  a  5° 
by  5°  grid  which  tends  to  average  the  fuie  structure  expected  in  the  cusp/cleft  area.  In  addition,  the 
observations  correspond  to  a  range  of  K.  indices,  while  the  modelled  results  are  calculated  for  Kp  =  3. 
Spndre  SttpmQord  is  at  the  location  of  the  dayside  boundary  between  the  polar  cap  and  the  sub  auroral 
region,  and  as  a  result,  neutral  wind  measurmnents  acquired  during  midday  periods  will  be  very  sensitive  to 
the  britiiHa  of  the  cusp/clefL  It  is  not  clear  from  this  data  set  the  cause  of  the  discrepancy  between  the  midday 
winds  Old  the  model,  but  it  seems  likely  that  Ae  differences  are  due  to  sampling  winds  in  or  about  the  cleft. 

In  sumnury,  we  have  sorted  several  years  of  optical  and  radar  data  based  upon  the  sign  of  the  IMF  By  com¬ 
ponent  In  all  cases,  the  IMF  B^  component  was  taken  to  be  near  zero,  or  (lefinitely  souAward.  The  binned 
averages  from  the  optical  measurements  m  Thule  and  Spndre  Sttpmfjord  have  been  displayed  m  a  polar  dial 
form  for  boA  signs  of  By  and  indicate  a  definite  asymmetry  m  the  cell  size  of  the  neutral  wind  convection 
pattern  for  By  positive  conditions.  The  cells  appearto  be  similar  in  sizeforthe  By  negative  case.These  obser- 
vationsareconsistentwiA  the  empirical  plasmaconvection  patterns  described  by //eppner  and  Afaynord/lO/. 
Neutral  horizontal  winds  m  the  F-region  obtained  hom  the  VSH  model  indicates  a  good  cotieq>ondence 
wiA  observations  acquired  at  Sbndre  StrpmQord  during  periods  m  which  this  station  is  wiAin  the  polar  cap. 
The  greatest  discrepancy  between  tpodeUed  and  measured  winds  occurs  during  the  cusp  transit  over  Sdndre 
Stromfjord. 
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at:  Space  Physics  Research  Laboratory,  Department  of  Atmospheric 
and  Oceanic  Science,  The  University  of  Michigan,  Ann  Arbor,  Ml 
48109-2143  USA) 

Measurements  of  F-region  neutral  winds  at  Sondre  Str0mQord  have 
been  routinely  made  by  a  Fabry-Perot  interferometer  since  early  1983. 
Incoherent  scatter  experiments  in  the  field  aligned  direction  have 
coincided  with  some  of  the  optical  observations.  In  particular,  neutral 
F  region  winds  have  been  extracted  from  76  radar  experiments 
totalling  2200  hours  observing  time  between  April  1983  and  July 
1988.  The  incoherent  scatter  data  analysis  has  shown  that  the  IMF  By 
component  has  a  significant  effect  on  the  meridional  wind:  for  By  < 
0,  equatorward  winds  are  enhanced  for  the  nightside  as  compared  to 
By  >  0  conditions.  The  current  study  supplements  the  previous 
analysis  by  extending  the  database  to  include  neutral  winds  measured 
by  the  Fabry-Perot  interferometer.  The  dependence  on  the  neutral 
meridional  wind  of  the  various  permutations  of  IMF  By  and  Bz 
positive  and  negative  will  be  presented  and  discussed  for  the  total 
database. 


Measurements  at  Thule.  Greenland 


Sequence  Using  Ground-Basec 


E.  S.  TRUDELL.  T.  L.  KILLEEN,  and  R.  J.  NICIEJEWSKI,  (All 
at:  Space  Physics  Research  Laboratory,  E>epartment  of  Atmospheric 
and  Oceanic  Science,  The  University  of  Michigan,  Ann  Arbor,  MI 
48109-2143  USA) 

SU.  BASU,  (Institute  for  Space  Research,  Boston  College,  Newton. 
MA  02159) 

B.  W.  REINISCH,  (University  of  Lowell  Center  for  Atmospheric 
Research,  450  Aiken  St.,  Lowell,  MA  01854) 

An  investigation  of  a  series  of  polar  cap  arcs  will  be  presentee  Data 
obtained  on  12  December  1988  during  a  High  Latitude  Plasma 
Structure  (HLPS)  campaign  will  be  used  to  discuss  the  structure  and 
dyriamic  behavior  of  the  upper  atmosphere  during  this  occurrence.  A 
series  of  OI  (6300  A)  and  OI  (5577  A)  all-sky  images  and  Fabry- 
Perot  interferometer  measurements  were  acquired  at  the  University  of 
Michigan  airglow  facility  in  Thule,  Greenland  (A=86®)  on  December 
12,  1988  from  0851  UT  to  1120  UT.  The  all-sky  images  depict 
intense  sun-aligned  arcs  later  follows  d  by  several  diffuse,  horseshoe¬ 
shaped  arcs.  Simultaneous  FPI  data  obtained  for  both  OI  emissions 
show  elevated  temperatures,  and  initially,  strong  anti-sunward  winds. 
Additional  data  has  been  obtained  through  digisonde  measurements  in 
the  form  of  electron  densities  and  ion  drifts,  and  irregularity  drift 
measurements  using  the  spaced  receiver  scintillation  technique.  These 
data  will  be  used  along  with  the  optical  signature  of  the  arcs  to 
de^ribe  the  morphology  and  dynamics  at  altitudes  of  1 10  Idlometers 
using  the  OI  (5577  A)  emission  and  220  kilometers  using  the  OI 
(6300  A)  emission. 
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Ground  Based  OI  (6300  A)  Fabrv  Perot  Interferometer  Observations 
from  Thuie  and  Sondre  Stromf^ord.  Greenland:  Svstematics  in  the  F 
Region  Neutral  Winds  Observed  Concurrently  bv  Both  Instruments 

R  J  NrriEJEWSKI.  J.  P.  THAYER,  and  T.  L.  KILLEEN  (AJl  at: 
Space  Physics  Research  Laboratory,  Department  of  Atmospheric  and 
Oceanic  Science,  The  University  of  Michigan,  Ann  Arbor,  MI  48109- 
2143  USA) 

Measurements  of  polar  upper  atmosphere  OI  (6300  k)  emissions  by 
Fabry  Perot  interferometers  at  both  Thule  and  Spndre  Strpmfjord, 
Greenland,  have  been  acquired  since  1986.  Previously,  studies  of  the 
systematics  of  F  region  neutral  winds  have  been  performed  using 
only  a  single  station  [eg.  Niciejewski  et  al,  EOS,  1989].  The  two 
stations  are  separated  by  ~9()0  km,  and  thus  during  clear  sky 
conditions,  both  Fabry  Perot  interferometers  can  simultaneously 
observe  a  large  portion  of  the  polar  cap  F  region.  In  this  presentation, 
we  will  discuss  results  from  two  modes  of  operation:  1)  concurrent 
observations  of  the  same  volume  from  both  stations  obtained  during 
the  1987/88  observing  season,  and  2)  independent  meridional  and 
zonal  neutral  wind  observations  obtained  at  the  same  time  from  both 
stations.  The  latter  will  be  compared  with  VSH  model  predictions  of 
the  neutral  winds  for  a  case  of  low  geomagnetic  activity. 


Ground-based  observations  of  ion/neutral  coupling  at  Thule  and 
Oanaq.  Greenland:  IMF  dependence 

L  E.  THAYER.  R.  J.  NICIEJEWSKI.  and  T.  L.  KILLEEN.  (All  at: 
Space  Physics  Research  Laboratory,  I^partment  of  Atmospheric  and 
Oceanic  Science,  The  University  of  Michigan,  Ann  Arbor,  MI  48 109- 
2143  USA) 

J.  BUCHAU,  (Geophysics  Lab,  Hanscom  AFB,  Bedford,  MA 
01731-5000) 

B.  W.  REINISCH,  and  G.  CROWLEY  (Both  at:  University  of 
Lowell  Center  for  Atmospheric  Research,  450  Aiken  St.,  Lowell,  MA 
01854) 

During  December,  1988,  twenty-four  hours  of  darkness  and  clear 
sky  conditions  permitted  continuous  observations  of  the  OI  (6300A) 
airglow  by  a  Fabry-Perot  interferometer  located  at  Thule,  Greenland. 
Thus,  a  long  term  continuous  record  of  the  F-region  neutral  winds 
was  obtained  for  that  month.  During  this  same  time  period,  a  digital 
ionosonde  located  at  Qan§q,  Greenland  (1 10  km  due  north  of  Thule) 
was  in  operation  measuring  electron  density  profiles  and  F-region  ion 
drifts.  Using  these  diagnostics,  a  one  to  one  correspondence  between 
the  F-region  neutral  wind  and  the  ion  drift  measurements  may  be 
established.  This  allows  for  the  investigation  of  ion/neutral  coupling 
from  the  ground-based  observations  at  a  temporal  resolution  of  about 
15  minutes.  From  the  lb**!  to  the  24th  of  December,  IMF  data  from 
the  IMP  8  satellite  were  available  and  indicated  intervals  of  B2  directed 
northward  IMF  conditions  during  this  period.  Signatures  of  sunward 
directed  neutral  winds  and  ion  drifts  were  observed  coincident  with 
northward  IMF.  The  response  of  the  neutral  winds  to  sunward  ion 
drifts  under  B2  northward  conditions  is  investigated. 
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Niciejewski,  R.  J.,  T.  L.  Killeen,  R.  M.  Johnson,  and  J.  Thayer,  The  behaviour  of  the 
high-latitude  F-region  neutral  thermosphere  in  relation  to  IMF  parameters 

Ground  based  incoherent  scatter  radar  and  Fabry  Perot  interferometer  studies  in  the 
northern  high  latitudes  during  the  period  1983  to  1989  have  shown  that  the  F-region  neutral 
wind  field  pattern  depends  upon  the  sign  of  the  IMF  parameters.  For  example,  the  cell 
structure  of  the  northern  hemisphere  high  latitude  neutral  wind  field  during  periods  of  low 
geomagnetic  activity  depends  to  a  large  degree  upon  the  sign  of  the  IMF  By  parameter. 
Long  term  monitoring  of  the  F-region  thermosphere  by  FPI  in  Thule,  Greenland,  and  by 
both  m  and  ISR  in  Spndre  StrpmQord,  Greenland,  have  made  it  possible  to  produce 
maps  of  average  meridional  and  zonal  wind  fields  for  various  IMF  configurations  for 
northern  high  latitudes.  Comparison  of  observations  with  theoretical  wind  field  modelling, 
such  as  the  Vector  Spherical  Harmonic  model,  indicates  that  most  observed  features  are 
within  the  context  of  the  models. 


Fabry  Perot  interferometer  observations  of  neutral  winds  in  the  polar 
cap  during  northward  IMF 

R.  J.  NICIEJEWSKI.  T.  L.  KILLEEN,  and  Y.  WON  (All  at:  Space 
Physics  Research  Laboratory,  Department  of  Atmospheric  and 
Oceanic  Science,  The  University  of  Michigan,  Ann  Arbor,  MI  48 109- 
2143  USA) 

Remote  observations  of  F-region  neutral  winds  have  been  performed 
by  the  University  of  Michigan  at  Thule  Air  Base  (lat:  76.5  °N,  long: 
69.0  ®W)  and  S0ndre  Str0mf|ord  (lat;  67.0  °N,  long:  50.9  “W), 
Greenland,  since  early  1983.  The  former  site  is  located  well  within 
the  polar  cap,  while  the  latter  is  situated  beneath  the  auroral  oval  dur¬ 
ing  dawn  and  dusk,  and  within  the  polar  cap  at  night.  Nearly  a  com¬ 
plete  solar  cycle  of  data  are  currently  available.  The  sign  of  the  Bj 
component  of  the  IMF  plays  an  important  role  in  the  dynamics  of  the 
polar  cap  ionosphere,  and  through  ion  drag  forcing,  the  polar  cap 
thermosphere.  Observations  indicate  that  the  anti-sunward  horizont^ 
neutral  wind  component  at  F  region  altitudes  is  much  reduced  during 
northward  IMF  conditions.  Examples  of  these  measurements  will  be 
shown  and  discussed  relative  to  standard  diuraally  reproducible  Vec¬ 
tor  Spherical  Harmonic  (VSH)  model  representations  of  the  neutral 
wind  field  for  250  km  altitude.  Averaged  horizontal  neutral  winds  for 
both  signs  of  B2  for  both  sites  will  also  be  displayed  and  discussed. 
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GROUND-BASED  OPTICAL  MEASUREMENTS  FROM  WITHIN  THE 
GEOMAGNETIC  POLAR  CAP  AT  THULE,  GREENLAND 

T.  L.  Killeen*.  R.  J.  Niciejewski,  Y.  Won, 
and  G.  S.  N.  Muny 
Space  Physics  Research  Laboratory 
Department  of  Atmospheric,  Oceanic,  and  Space  Sciences 
The  University  of  Michigan 
Ann  Arbor,  Michigan  4S109-2143 


Observations  of  the  dynamics  and  thermodynamics  of  the  polar  cap  thermosphere 
have  been  obtained  routinely  since  the  winter  of  1984  with  an  automated  optical 
observatory  located  at  Thule,  Greenland  (A®86®).  The  instrumentation  at  the 
observatory  includes  a  Fabry-Perot  interferometer  (FT^I),  a  half-meter  Ebert-Fastie 
spectrophotometer  (EFS),  and  a  digital  all-sky  camera  (ASC)  system.  The  FPl  has 
a  10  cm  etalon  coated  for  operation  between  SS77-A  and  7320-A  and  an  image 
plane  detector  similar  to  that  flown  on  the  Dynamics  Explorer  spacecraft  The  FPI 
and  EFS  instruments  are  controlled  automatically  throughout  each  8-mcmth  winter 
observing  season  (September  -  April)  with  an  LSI  11/23  computer  system  and 
CAMAC  interface  cards.  The  ACS  is  operated  on  a  campaign  basis.  The 
thermospheric  neutral  winds  and  temperatures  measured  with  the  FPI  in  the  years 
1985- IS^l  show  strong  dependencies  on  the  solar  cycle,  local  time,  geomagnetic 
activity  level  and  interplanetary  magnetic  field  orientation.  The  E^  has  been  u^ 
to  observe  various  nightglow  emissions  including  NI  (5200-A),  OI(5577-A), 
OI(6300-A),  Balmer  Ha(6563-A),  On(7320-30-A),  and  the  OH  (8-3)  vibrational 
band.  The  ACS  has  been  used  in  conjunction  with  the  FPI  to  correlate  sun-aligned 
polar  cap  arcs  and  thermospheric  d^amics.  This  paper  will  describe  the  current 
status  and  plans  for  the  automated  optical  observatory  at  Thule  and  will  discuss 
several  of  the  more  significant  scientific  findings  made  over  the  past  several  years. 
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OBSERVATIONS  OF  NEUTRAL  WINDS  IN  THE  POLAR  CAP 
DURING  NORTHWARD  IMF 

by 


R.  J.  Niciejewski.  T.  L.  Killeen,  and  Y.  Won 
Space  Physics  Research  Laboratory, 

The  University  of  Michigan 
Ann  Arbor,  Ml,  USA  48 109 

Long  term  remote  observations  of  neutral  winds  at  F-region  altitudes  have  been  per¬ 
formed  at  Thule  Air  Base  (lat:  76.5®  N,  long:  69.0  W),  Greenland,  and  Spndre  StrpmijQoid 
(lat:  67.0  N,  long:  50.9  W),  Greenland,  since  1983.  The  former  site  is  very  close  to  the 
geomagnetic  pole,  while  the  latter  site  is  within  the  polar  cap  for  several  hours  each  night 
on  either  side  of  geomagnetic  midnight.  The  data  base  that  currently  exists  corresponds  to 
nearly  one  full  solar  cycle.  This  paper  will  present  the  following:  one  hour  averages  of 
thermospheric  neutral  winds  for  northward  IMF  conditions  obtained  by  both  instruments 
binned  according  to  the  sign  of  By;  individual  data  sets  from  Thule  showing  the  depen¬ 
dence  of  the  neutral  wind  on  the  sign  of  B2;  and  deviations  of  the  averaged  northward  IMF 
measurements  from  standard  diuraally  teproducible  Vector  Spherical  Harmonic  (VSH)  rep¬ 
resentations  of  the  polar  cap  neutral  wind  field  for  250  km  altitude.  Thermospheric  wind 
measurements  at  Thule  Air  Base  during  northward  IMF  very  often  show  a  dramatic  reduc¬ 
tion  in  the  anti-sunward  horizontal  wind  component.  Examples  of  these  measurements  will 
be  shown. 
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and  Winds  During  the  1991/92  Optical  Observing  Season 

R.  J.  NICIEJEWSKI,  T.  L.  KILLEEN.  Y.  WON.  and  M.  TURN- 
BULL  (All  at:  Space  Physics  Research  Laboratory.  Department  of 
Atmospheric  and  Oceanic  Science,  The  University  of  Michigan,  Ann 
Arbor,  MI  48109-2143  USA) 

Fabry-Perot  interferometer  observations  of  the  thermosphere  are  cur¬ 
rently  being  performed  from  several  locations  in  the  northern  hemi¬ 
sphere  by  an  aeronomy  group  at  the  University  of  Michigan  in  sup¬ 
port  of  several  efforts:  the  CEDAR  program,  the  UARS  ground  sup¬ 
port  program,  and  the  SPIRIT  II  ground  support  effort.  This  current 
observing  season,  late  August  '91  to  late  April  '92,  coincides  with  an 
apparent  increase  in  solar  activity,  and  the  subsequent  heating  of  the 
terrestrial  thermosphere  and  increase  in  geomagnetic  activity.  Neutral 
winds  (both  vertical  and  horizontal)  and  neutral  temperatures  are  rou¬ 
tinely  obtained  from  these  sites,  which  encompass  a  wide  range  of 
latitudes  and  local  times.  In  this  presentation,  the  current  set  of  geo¬ 
physical  measurements  will  be  compared  with  previous  measurements 
from  the  decline  phase  of  the  previous  solar  cycle  (1983  onwards), 
and  with  the  current  state-of-the-art  VSH  model.  In  addition,  the  util¬ 
ity  of  combining  multi-station  observations  on  single  "weather  maps " 
will  be  presented  for  various  events,  including  the  2/8  -  2/9,  1992, 
geomagnetic  storm  (Ap  =  65,  Kp  <  7). 


Seasonal  and  Diurnal  Variability  in  the  Northern  High  Latitude  Ther¬ 
mosphere:  Neutral  Temperatures  and  Horizontal  Winds 

R.  J.  NICIEJEWSKI.  T.  L.  KILLEEN,  and  Y.  WON  (All  at:  Space 
Physics  Research  Laboratory,  Department  of  Atmospheric  and 
Oceanic  Science.  The  University  of  Michigan,  Ann  Arbor,  MI  48109- 
2143  USA) 

Ground  based  measurements  of  the  neutral  wind  and  neutral  tempera¬ 
ture  in  the  thermosphere  have  been  made  for  both  solar  minimum  and 
solar  maximum  conditions  from  two  stations  in  the  high  northern  lati¬ 
tudes:  Thule  Air  Base  and  Sendre  StrpmQord,  Greenland.  For  this 
study,  which  is  a  part  of  the  Coordinated  Analysis  of  the  Thermo¬ 
sphere  effort  (CAT),  the  data  sets  have  been  screened  to  include  only 
periods  of  low  geomagnetic  activity  (Kp  <  3)  and  have  been  hourly 
binned  for  two  month  intervals  centered  on  the  autumnal  equinox,  the 
winter  solstice,  and  the  vernal  equinox.  Both  the  averaged  winds  and 
the  averaged  temperatures  show  interesting  seasonal  and  diurnal 
variations.  The  thermosphere  above  Sondre  StromQord  is  cooler 
during  the  winter  season  compared  to  the  equinoxes,  the  meridional 
winds  are  less  in  magnitude  during  the  winter,  and  the  zonal  winds 
are  similar  for  all  seasons.  Variability  in  the  seasonal  behavior  above 
Thule  Air  Base  is  not  as  marked  in  contrast  to  that  above  Sdndre 
StromQord.  However,  the  diurnal  variability  in  the  thermosphere 
above  Thule  Air  Base  is  more  pronounced  than  that  above  Sondre 
StrpmQord:  both  the  zonal  and  ^e  meridional  winds  exhibit  diurnal 
variations  nearly  double  for  the  northern  location  compared  to  the 
southern  site.  The  binned  data  will  be  presented  with  least  square  fits 
describing  both  the  diurnal  and  the  semi-diurnal  components,  as  well 
as  compared  to  recent  numerical  models,  and  the  underlying  driving 
processes  will  be  discussed. 
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Fabry-Perot  Interferometric  Measurements  of  Mesopause  and  Lower 
Thermospheric  Structure 

T.  L.  Killeen  and  R.  J.  Niciejewski  (Space  Physics  Research 
Laboratory,  The  University  of  Michigan,  Ann  Arbor,  MI  48109; 
313-747-3430) 

Fabry-Perot  interferometric  techniques  have  been  applied  in  recent 
years  to  investigate  the  thermal  and  dynamical  structure  of  the 
mesopause  and  lower  thermosphere  from  the  ground  as  well  as  from 
space.  This  paper  provides  a  review  of  the  new  instrumental 
approaches  and  geophysical  results  gained  using  the  Fabry-Perot 
technique.  Comparisons  with  numerical  and  empirical  models  will  be 
shown.  Emphasis  will  be  placed  on  the  results  obtained  with  the 
University  of  Michigan's  chain  of  opdeal  observatories. 
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Observations  of  high-latitude  thermospheric  neutral  winds  over  a  half 
solar  cycle  from  Thule  Air  Base  and  Sandre  Strflmfiord.  Greenland 

X.  WQM.  T.  L.  KILLEEN,  and  R.  J.  NICIEJEWSKI  (All  at:  Space 
Physics  Research  Laboratory,  Department  of  Atmospheric  and 
Oceanic  Science,  The  University  of  ^chigan,  Ann  Arbor,  MI  48 109- 
2143  USA) 

Ground-based  Fabry-Perot  interferometers  (FPIs)  at  Thule  Air  Base 
and  Spndrestrpm,  Greenland  have  been  used  to  determine  the 
dynamical  and  thermodynamical  properties  of  the  thermosphere  by 
measuring  emission  line  [Hoflles  from  the  0( ^D)  and  0(^S)  nightglow 
at  high  spectral  resolution.  Data  has  been  obtained  routinely  since  the 
winter  of  1984  with  the  FPIs  at  two  automated  optical  ob^rvatories 
located  at  Thule  Air  Base  and  Spndie  StromQord,  Greenland.  The 
multi-year  data  base  has  been  analyzed  to  determine  the  average 
depen^nce  of  the  high-latitude  thermospheric  circulation  pattern  as  a 
function  of  solar  cylcle,  Kp,  IMF  orientation,  and  solar  F10.7.  The 
observed  winds,  obtained  over  a  period  corresponding  to  half  a  solar 
cycle  (1985  -  1^1),  have  been  compared  with  the  predictions  of  the 
VSH  and  HWM  models  of  thermosphere  dynamics.  This  period 
encompasses  both  solar  minimum  and  solar  maximum  conditions.  We 
present  the  averaged  data  sets  from  the  combined  Thule  and 
S0ndrestr0m  observations  and  discuss  them  in  the  context  of  the 
model  predictions. 
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Studies  of  High-Latitude  Lower  Thermospheric  Dynamics  Based  on 
Measurements  and  Modelling  of  Lower  Thermospheric  Coupling 
Study  Experiments 

R.  M.  Johnson.  Y.  Won,  R.  Niciejewski,  and  T.  L.  Killeen  (Space 
Physics  Research  Laboratory,  The  University  of  Michigan,  Ann 
Arbor.  MI  48109;  313-747-3430) 

T.  S.  Viidi  and  P.  J.  S.  Williams  (both  at  Physics  Dept,  University 
College  Wales,  UK) 

A.  H.  Manson  and  C.  E.  Meek  (both  at  ISAS,  University  of 
Saskatchewan,  Canada) 

R.  A.  Vincent  (University  of  Adelaide,  Australia) 

G.  H.  Fraser  (University  of  Canterbury.  New  Zealand) 

J.  Forbes  (Dept,  of  Electrical,  Computer,  and  Systems  Engineering, 
Boston  University,  MA) 

C.  Fesen  (Thayer  School  of  Engineering,  Dartmouth  College, 
Hanover,  NH) 

R.  Roble  (High  Altitude  Observatory,  NCAR,  Boulder,  CO) 

Observations  obtained  by  high-latitude  radars  and  optical  instruments 
during  Lower  Thermospheric  Coupling  Study  experiments,  in 
conjunction  with  modelling  studies,  allow  us  to  make  significant 
advances  in  our  understanding  of  the  structure  and  dynamics  of  the 
lower  thermosphere.  Seasonal  variations  in  high-latitude  lower 
thermospheric  mean  winds  and  temperatures  and  diurnal  and 
semidiurnal  "tidal"  components  are  apparent  in  observations.  During 
periods  of  intense  geomagnetic  activity,  evidence  for  accelerated 
neutral  flows  in  the  lower  thermosphere  can-  be  found.  Comparison 
of  observational  results  with  TIGCM  simulations  of  LTCS 
experiments  generally  shows  some  agreement  with  the  features  seen 
in  the  data,  although  the  magnitude  of  the  features  are  generally 
smaller  and  the  timing  slightly  different  in  the  model  results. 
Semidiurnal  tidal  amplitudes  determined  from  high-latitude  LTCS 
measurements  are  generally  larger  than  those  calculated  using  tidal 
models,  and  there  are  significant  variations  in  observed  and  mc^elled 
phase  profiles. 
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Seasonal  variation  of  winds  in  the  high  latitude 
lower  thermosphere/upper  mesosphere  at  Thule 
Air  Base  and  Sondre  Stromfjord,  Greenland. 

Y  Won.  T  L  Killeen,  R  M  Johnson  and  R  J  Niciejewski 
(Space  Physics  Research  Laboratory,  Department  of 
Atmospheric,  Oceanic,  and  Space  Sciences,  The 
University  of  Michig.in,  Ann  Arbor,  MI  48109-2143) 

Ground-based  Fabry-Perot  interferometers  (FPI) 
have  been  used  to  determine  dynamical  properties  of  the 
lower  thermosphere/upper  mesosphere  by  measuring 
emission  line  profiles  from  the  0(^S)  nightglow  at  high 
spectral  resolution.  The  resulting  winds  are  normally 
referred  to  the  altitude  of  the  peak  of  the  0(^S)  emission 
layer,  near  97  km.  Data  has  been  obtained  since  the  winter  of 
1987  at  two  optical  observatories  located  at  Thule  Air  Base 
and  Sondre  Stromfjord,  Greenland.  In  our  analysis,  data 
have  been  removed  during  periods  of  cloud  cover  and  high 
signal  intensity  to  eliminate  auroral  contaminations.  Average 
neutral  winds  are  calculated  and  spectral  analysis  is  cairi^ 
out  to  examine  the  major  tidal  wind  components.  The 
seasonal  variations  of  these  wind  components  will  be 
discussed  and  compared  to  model  predictions. 
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Y.  WON.  R.  M.  JOHNSON.  T.  L.  KILLEEN,  and  R.  J. 
NlCIEJEWSKl  (All  at:  Space  Physics  Research  Laboratory. 
Department  of  Atmospheric  and  Oceanic  Science,  The  University  of 
Michigan,  Ann  Arbor,  MI  48109-2143  USA) 

J.  M.  FORBES  (High  Altitude  Observatory.  National  Center  for 
Atmospheric  Research,  Boulder,  CO,  80307) 

Measurements  of  winds  in  the  lower  thermopshere  and  upper 
mesosphere  region  were  obtained  from  various  ground-based  Fabry- 
Perot  interferomteters  and  incoherenet  scatter  radar  observatories. 
Tidal  information  was  then  extracted  by  performing  harmonic  analysis 
on  these  data.  The  day-to-day  variation  of  these  tidal  fields  are 
discussed  in  this  presentation,  in  addition  to  our  description  of  the 
seasonal  characteristics.  The  amplitudes  and  phases  of  the  seasonally 
averaged  data  are  utilized  to  obtain  a  empirical  description  of  the 
semidiurnal  tidal  fields  by  using  a  set  of  numerically  generated 
"Hough  Mode  Extension"  (HI^)  functions  for  u,  v,  and  St 
corresponding  to  the  (2,2),  (2,3),  (2,4),  and  (2,5)  semidiurnal  tidal 
modes.  The  results  will  he  compared  with  tidal  fields  from  the  Forbes 
and  Vial  (1989)  numerical  tidal  model. 

Thermflaphfidc  Neutral  Wind  Climatology  in  the  Northern  High 

Latimdss 

R.  J.  NlCIEJEWSKl  (Space  Physics  Research  Laboratory, 
Department  of  Atmospheric  and  Oceanic  Science,  The  University  of 
Michigan,  Ann  Arbor,  MI  48 109-2 143  USA) 

J.  P.  THAYER  (SRI  International,  Geoscience  and  Engineering 
Center,  Menlo  Park,  CA  94025) 

A.  L.  Aruliah  (University  College  London,  England) 

M  part  of  the  Coordinated  Analysis  of  the  Thermosphere  (CAT),  a 
high-latitude  working  group  was  formed  to  analyze  and  interpret  the 
diumal/seasonal/solar  cycle  variations  in  the  observed  thermospheric 
neutral  wind  during  geomagnetically  quiet  conditions.  The  high- 
latitudes  sites  which  have  partiepated  in  the  study  include  Thule, 
Greenland  (76.5N,  -68.4E),  Sondrestrom,  Greenland  (67.0N,  - 
50.9E),  and  Kiruna,  Sweden  (67.8N,  20.4E).  The  Thule, 
Sondrestrom,  and  Kiruna  sites  provide  full  vector  winds  from  Fabry- 
Perot  interferometer  (FPI)  measurements  of  the  thermospheric  6300A 
emission  while  thermospheric  neutral  winds  in  the  magnetic  meridian 
are  derived  from  the  incoherent-scatter  radar  (ISR)  measurements  at 
Sondrestrom.  A  large  database  of  thermospheric  neutral  winds  has 
been  complied  covering  an  entire  solar  cycle  (1981  -  1991).  TTie  ITI 
data  are  limited  to  nighttime  observations  and  therefore  exclude  the 
summer  season  as  well  as  some  daytiiiK  hours  throughout  the  year. 
Comparisons  of  these  data  sets  with  the  NCAR  TIGCM  simulation 
and  the  HWM90  empirical  model  will  be  presented.  A  harmonic 
analysis  was  applied  to  the  averaged  ISR  meridional  winds  to 
demonstrate  the  daily  variability  and  the  seasonal  dependence  of  the 
fitted  harmonics.  Care  must  be  taken  in  interpreting  the  harmonic 
ai^ysis  as  tidal  in  nature  as  the  forcing  by  tlw  two-cell  ion  convection 
will  introduce  a  semidiurnal  influence  on  the  observations. 
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NEUTRAL  WIND  VORTICES  IN  THE  HIGH-LATITUDE 
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